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ABSTRACT 


Nonlinear regression analysis holds significant popularity in mathematical, 
engineering, and social science domains. Disciplines like financial matters, 
biology, and natural chemistry have broadly utilized nonlinear regression 
models (NLRMs). Cloned datasets have their own importance in such areas 
which provide the same fit of bivariate and multivariate nonlinear 
regression models for the actual datasets. This article presents a sequence 
of cloned datasets that give exactly the same fit of bivariate and multivariate 
nonlinear regression models. 


Keywords: cloned data, nonlinear regression model, fictitious datasets, 
data visualization. 


1. INTRODUCTION 


If genuine information is private and cannot be shown, a matching or 
alternative set of data is required which provide same summary Statistics as 
of the actual data. Cloned data refers to the alternative or matching set of 
data through mathematical techniques that allow rapid provisioning in 
testing and developments. Data cloning has its own significance as an 
alternative method for protecting confidential information and database. 
Table 1 shows four fictitious distinct cloned datasets (CDSs) created by 
Anscombe [1] to demonstrate the significance of graphs in statistical 
analysis. The summary statistics (mean, standard deviation, and correlation) 
as well as the parameter estimates of the fitted regression equation R? and 
estimated standard deviation of residuals are identical across these four 
distinct CDSs, however, they were vastly different scatter plots as shown in 
Figure 1. Dataset I was strongly linear with a single outlier and II appears 
to follow a parabolic distribution, whereas dataset III appears to adhere to 
a noisy linear regression model (LRM), and dataset IV appears to follow a 
vertical line with the regression thrown off by a single outlier. Datesets in 
Table 1 are significant and frequently used to show how important visible 
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methods are. These datasets were also known for their significant use in 
education. However, the method used to create the datasets was not 
explained in [1]. A genetic algorithm-based approach was proposed by 
Chatterjee and Firat [2], who generated 1,000 random datasets with 
comparable summary statistics and graphics for the basic datasets. 
Govindaraju and Haslett [3] devised a method for producing datasets by 
regressing the response on the covariate in the direction of their 
unconditional sample means, while maintaining identical LRM estimates. 
As a result, the variability in the response and the covariate decreased in 
each subsequent cloned dataset. Haslett and Govindaraju's [4] method for 
creating matched datasets was extended to include a multiple linear 
regression model, ensuring that the matched datasets have an identical fit to 
the original data. The idea of data-cloning emerged from both biostatistics 
[5, 6] and financial time series [7]. 


Cloning for maximum likelihood estimation using Bayesian software 
was achieved by the simple device of replicating the original data many 
times [6]. Fung et al. [8] expressed that the creation of CDSs to anonymize 
sensitive data was another application for datasets with the same statistical 
properties, as discussed in [3]. In this instance, it is critical that individual 
data points were altered, while the data's overall structure remained 
unchanged. 


Haslett and Govindaraju [9] described a straightforward approach for 
modifying LRM data, while still obtaining the same fitted regression 
parameters. Ponciano et al. [10] showed how structural parameter non- 
identifiability can be diagnosed with Data Cloning (DC) and distinguished 
from other parameter estimability issues, such as when parameters are 
structurally identifiable but not estimable in a given data set or when they 
are identifiable and weakly estimable. Bayesian phylogenetics software can 
be used to diagnose non-identifiability with the DC approach. Additionally, 
it was demonstrated that DC can be used to examine and eliminate the 
influence of priors, particularly when prior elicitation was difficult. Finally, 
DC can be used to investigate at least two significant statistical issues when 
applied to phylogenetic inference, developing effective sampling strategies 
for computationally expensive posterior densities, and evaluating the 
identifiability of discrete parameters, such as the tree's topology. 


Data confidentiality is one of the designed goals of tunable encrypted 
deduplication, see Amvrosiadis and Bhadkamkar [11]. Additionally, it 
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reduced the risk of data leakage brought by frequency analysis. 
Furthermore, it was identified that better ways of seeing and exploring data 
lead to better insights. The "Datasaurus" Cairo dataset was created by 
Alberto Cairo [12]. This, like Anscombe's Quartet, emphasized the 
significance of data visualization, despite the dataset's normal summary 
statistics, the plot it produced depicted a dinosaur. They started with the 
datasaurus and created additional datasets with the same summary statistics. 
Additionally, Cairo's Datasaurus data visualization prohibited to solely rely 
on the summary statistics of the used data. 


Resultantly, according to [2], datasets should be as graphically distinct 
as possible. With different standard deviations but identical means and 
LRM estimates, [3, 4, 9] data are intended to be graphically comparable. 
Matejka and Fitzmaurice [13] developed a novel method for creating 
datasets, which are identical across a variety of statistical properties but 
visually distinct during the data exploration. To address the primary 
empirical facts of financial time series, numerous complex parametric 
stochastic volatility models were proposed in the subsequent literature. The 
models that Mao et al. [14, 15] proposed incorporated a broader 
asymmetric volatility function. 


Hussain et al. [16] used a simple procedure to clone data for nonlinear 
regression models having linearizable or nonlinearize regression functions, 


such as aX?, ab, ae, kab”, ks* bo", ktabX, —“_, A [axy? +(1- 


1+bcX’ 


1 
a)X;°] ?. They found that cloned data generated by linearizable or non- 
linearizable estimable functions of parameters have unchanged estimates. 
The procedure increased the sample size of cloned data without changing 
the parameters estimates, which was for n original sample points (x, y). This 
generated n? observations by adding [aj: i= 1, 2, ..., n] to the data points y 
over );a; = 0. Due to increased sample size, cloned estimates showed 
smaller standard errors as compared to the original standard errors. This 
procedure used by [16] was sufficient for the first iteration because in the 
next iterations, it became tedious work. This procedure was useful for 
modeling but not for confidentialising or encrypting data, as in the design 
matrix variables remained unchanged. In this case, the term 
"confidentializing" referred to making the values of particular variables 
certain for particular people that cannot be deduced from the data. Our goal 
in this article is to create datasets with the same fit for nonlinear linear 
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regression models (NLRMs). [3, 4] methods were used to generate these 
cloned data sets. To get around the problem in [16], nonlinear regression 
models with linearizeable regression functions were the prime focus of this 


article. 


Table 1. Anscombe’s CDSs with Pairs (x, y1), (x, y2), (x, y3) and (x4, y4) 


x yı y2 y3 x4 y4 

7 6.42 7.26 4.82 8 7.91 
4 5.39 3.10 4.26 19 12.5 
11 7.81 9.26 8.33 8 8.47 
13 12.74 8.74 7.58 8 7.71 
8 6.77 8.14 6.95 8 5.76 
5 5.73 4.74 5.68 8 6.89 
12 8.15 9.13 10.84 8 5.56 
6 6.08 6.13 7.24 8 5.25 
10 7.46 9.14 8.04 8 6.58 
9 7.11 8.77 8.81 8 8.84 
14 8.84 8.10 9.96 8 7.04 


Fitted Line Plot 
Y=3.002+0.4997X 


Fitted Line Pilot 
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Figure 1. Scatter Plots of Anscombe’s Datasets with Matched Simple 


Regression Models 
2. THE NONLINEAR REGRESSION MODEL 
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2.1. The Regression Model (RM) 
The Regression Model (RM) talks about the relationship between a 
response variable Y and one or more covariates XV). The general model is: 


Y, = HOG XA, ai Aia seip) +E; 


Here, H is a suitable function that relies on the covariates 
LIR aa and parameters A1,A2, .., Ap . The unstructured 
deviations from H are defined by means of random errors (REs) 
E;~N(0,07). 


2.2. The Linear Regression Model 


In multiple LRM (MLRM), function H are characterized as linear in the 
parameters. 


sie ae wun XO: Ay, Ay, aA) SAX $A, x Pacts” 


where X 7 ) can be arbitrary functions of the original covariates X ~ ), 
2.3. The Nonlinear Regression Model 


In NLRM, function H is regarded in such a way that it can’t be written 
as linear in parameters. In case, there are infinite ways to explain the 
deterministic part of the model. 


2.4. Linearizable Regression Functions (LRFs) 


In NLRMs, functions Hcan be linearized by the transformation of the 
variable of interest and the explanatory variables. Therefore, the regression 
is named as function H which is linearizable if it can be converted into a 
function linear in the parameters. 


2.5. A Few Examples of Nonlinear Regression Functions 
1- H(X;; A41, A42) = aa” 

2- H(X; A1, A2) = A, Ax 

3- H(X; Ay, A2) = A,e42*: 

4- H(X®, XP; Ay, Az As) = Ay (XP) 42 (x) Aa 


2.6. Linearizable Regression Function Model 
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A LRM with the LRF in the referred example is based on the model 
given below: 


In(¥) = By + BX? + BX + By XO = In(x), XP = n(x) 


Where, E; follows the normal distribution. This model was back-converted 
and for this reason, the following equation was obtained: 


Y; = Ay (X)42(x) 438, ; B, = exp(E;),i = 1,2,....1 


The errors E; follows lognormal distributed and contributed 
multiplicatively. The assumptions about the random deviations were 
accordingly now appreciably distinct for a model, which was primarily 
based on: 


Y; = A (Xp?) 2X) 49 + BF 

with random deviations E; that follows normal distribution and contributed 
additively. 

3. DATA CLONING BY USING REGRESSING Y ON X AND X ON 
Y 


Assuming n paired observations of X and Y say (x;, yi) i = 1,2,°::,n. 
The following procedure from [3] would generate a sequence of CDSs by 
obtaining the same fitted NLRM equations. 


3.1. Procedure for Bivariate Nonlinear Regression Model Y = AX® 


The simple NLRM (a geometric or power curve) Y = AX? was 
linearizable due to logarithmic transformation as Y = a + BX where Y = 
In(Y),X = In(X),a = In (A) ,and A = exp (a). The inverse nonlinear 

1 
regression model (INLRM) of Y = AX? is = (=P , which was also 


mO p 
D = 


linearizable as X=c+dY where d= =,D = 2 ,c=— 
È 

exp (— 7). 

1. First fit regression Y on X, namely Y, = a + B. Also fit inverse 
regression (IR) Ž on Y, namely X, = c + dľ. 

2. The regression of Y, on X, would be Y> = a + BX,, preserving same 
parameter estimates. Likewise, X, = c + dY,. Note that S$ < S$ and 
Sg, < SZ. 
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3. The above method can be iterated with Y, and Š, as done in step 1 to 
gain cloning sets of data having the identical linearizable regression 
equation (LRE). Again S$, < SS < Sg and Ss, < Sg, < SZ and so on. 

4. If preferred, transform back for CDSs, having same coefficients of 
NLRM. It was noted that variability in Y and X of the cloned datasets 
fluctuated after every iteration (see Table 2). 


Example 1. Consider the variables X= (0.5, 1.5, 2.5, 5.0, 10.0)' and Y= 
(3.4, 7.0, 12.8, 29.8, 68.2)", resulting in the nonlinear regression fit 


f = 5.709057X1-01876 


Steps 1-4 given above would yield the CDSs proven in Table 2. with exactly 
the same equation of fitted NLRM as in (Eq. 3.1). 


(3.1) 


Table 2. Cloned Data Sets Having the same Non Linear Regression Fit Y = 


AX? 
Raw data First iteration Second iteration 
X Y Xi Yı X2 Yo 
0.5 3.4 0.61884 2.81765 0.51656 3.50135 
1.5 7.0 1.23925 8.62897 1.51542 7.10350 
2.5 12.8 2.21417 14.52011 2.49958 12.83071 
5.0 29.8 4.99031 29.42030 4.92915 29.36222 
10.0 68.2 11.06349 59.61073 9.72023 66.07545 
Mean 3.9 24.24 4.02521 22.99955 3.83619 23.77465 
Variances 14.425 706.448 18.2784 516.8321 13.50223 657.32050 
Correlation - 0.99684 - 0.99803 - 0.99679 
Third iteration Fourth iteration Fifth iteration 
X3 Y3 X4 Y4 X5 Y5 
0.63657 2.91275 0.53332 3.60356 0.65443 3.00903 
1.25687 8.71936 1.53068 7.20638 1.27437 8.80883 
2.21927 14.51765 2.49918 12.86087 2.22429 14.51524 
4.91979 28.99566 4.86072 28.93958 4.85168 28.58561 
10.73187 57.91217 9.45374 64.05832 10.41665 56.29510 
Mean 3.952877 22.61152 3.77553 23.33374 3.88428 22.24276 
Variances 17.04036 483.35930 12.65003 612.24240 15.9017 452.4788 
Correlation - 0.99797 - 0.99672 - 0.99791 


3.2. Procedure for Bivariate Nonlinear Regression Model Y = AB* 


A simple nonlinear regression model (an exponential curve) Y = AB* 
was linearizable due to logarithmic transformation as Y = a + bXwhere 


s —S 


mmi 
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Y = In(Y),a = In(A),A = exp (a), b = In(B), and B = exp (b). The 
1 
inverse nonlinear regression model of Y = AB* was X = In()*o,, which 


was also linearizable, as X = c + dY where d = 


In (C) = _ ©. 
In (D) „C = exp ( D 


1 1 
noy? = exp -c = 


1. First fit the regression of Y on Xspecifically Ý, = a + bX. Also, fit the 
IR of X on Ÿ particularly X, = c + dř. 

2. The regression of Y, on X, would be Y, =a+bX,, maintaing the 
identical parameter estimates. Similarly, X = c + dY,. Note that S$, < 
S$ and SZ, < Sf. 

3. The above prcedure can be iterated with Y, and X, as in step 1 to obtain 
CDSs having the same LRE. Again SẸ, < SẸ, < S and S$, < S$, < S 
and so on. 

4. If preferred, convert back to produce a sequence of CDSs, all with the 
same NLRM coefficients. Therefore, it was observed that variability in 
Y of the cloned datasets fluctuated after each and every generation (see 
Table 3). 


Example 2. Consider the variables X= (0, 1, 2, 3, 4, 5, 6, 7, 8)" and Y= 
(0.75, 1.20, 1.75, 2.50, 3.45, 4.70, 6.20, 8.25, 11.50)" and resulting in the 
nonlinear regression fit 


Ê = (0.8573324)(1.392474)* (3.2) 


Steps 1-4 described above would generate the CDSs presented by Table 3 
having exactly same NLRM fitted equation as in (Eq. 3.2). 
Table 3. Cloned Data Sets Having the same Nonlinear Regression Fit Y = 


AB* 
Raw data First iteration Second iteration 
X Y Xi Yı X2 Yo 
0 0.75 -0.3781 0.8573 0.0235 0.7565 
1 1.20 1.0332 1.1938 1.0176 1.2070 
2 1.75 2.1660 1.6624 2.0118 1.7563 
3 2.50 3.2370 2.3148 3.0059 2.5037 
4 3.45 4.2041 3.2233 4.0000 3.4486 
5 4.70 5.1325 4.4883 4.9941 4.6896 
6 6.20 5.9642 6.2499 5.9882 6.1762 
School of Science im. 
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Raw data First iteration Second iteration 

7 8.25 6.8219 8.7028 6.9824 8.2045 

8 11.50 7.8192 12.1184 7.9765 11.4143 
Mean 4 4.47778 4 4.53456 4 4.46187 
Variances 7.5 12.95069 7.45591 14.67662 7.41209 12.73009 
Correlation - 0.95442 - 0.91709 - 0.95498 

Third iteration Fourth iteration Fifth iteration 

X3 Y3 X4 Y4 X5 Ys 


-0.3524 0.8640 0.0469 0.7629 -0.3268 0.8707 
1.0506 1.2008 1.0352 1.2140 1.0679 1.2078 
2.1768 1.6688 2.0234 1.7626 2.1875 1.6753 
3.2415 2.3193 3.0117 2.5075 3.2459 2.3238 
4.2029 3.2233 4.0000 3.4473 4.2017 3.2233 
5.1258 4.4796 4.9883 4.6793 5.1192 4.4709 
5.9526 6.2256 5.9766 6.1526 5.9412 6.2016 
6.8053 8.6521 6.9648 8.1596 6.7888 8.6021 
7.7968 12.0245 7.9531 11.3298 7.7744 11.9318 


Mean 4 4.51756 4 4.44617 4 4.50081 
Variances 7.36852 14.41740 7.32521 12.51401 7.28215 14.16366 
Correlation - 0.91777 - 0.95553 - 0.91844 


3.3. Procedure for Bivariate Nonlinear Regression Model Y = Ae®* 


The simple nonlinear regression model (an exponential curve) Y = 
Ae®* was linearizable due to logarithmic transformation as Y = a + BX, 
where Y = In(Y),a = In(A), and A = exp (a). The inverse nonlinear 


1 
regression model of Y = Ae®* was X = nÓ , which was also 


linearizable as X=c+dY where denpa gasa 
E 
exp (— 7). 


1. First fit simple LRM of Y on X representing Y, = a + BX. Also, fit the 
simple inverse LRM of X on Y describing asX, = c + dY. 

2. The regression of Y, on X, would be Y, = a + BX;, saving the alike 
parameter estimates. In the same way, X> = c + dY,. Note that S$ < 
S$ and SZ, < SZ. 

3. The approach used above can be iterated with Y, and X, as in step 1 to 
get CDSs with the same LRE. Again SẸ, < S$, < S$ and S$, < SZ, < 


S¥ and so on. 
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4. If preferred, transform back to get a sequence of CDSs, all with the same 
NLRM coefficients. It was observed that variability in Y of the cloned 


datasets fluctuate after every iteration (see Table 4). 


Example 3. Consider the variables X= (0.5, 0.8, 1.4, 2.0, 2.5)" and Y= (9.1, 
8.5, 7.5, 6.7, 6.1)! and resulting in the nonlinear regression fit 


Ê = 9.989049e~0.1991399% 


Here, the CDSs would be yielded as shown in Table 4, using steps 1-4 
discussed earlier, to produce same equation of fitted NLRM as given in 


(3.3) 


(Œq. 3.3). 
Table 4. Cloned Data Sets Having the Same Nonlinear Regression Fit Y = 
AeB* 
Raw data First iteration Second iteration 
X Y Xı Yı X2 Y2 
0.5 9.1 0.46924 9.04235 0.50111 9.09792 
0.8 8.5 0.81135 8.51797 0.80076 8.49874 
1.4 75 1.43912 7.55866 1.40005 7.50000 
2.0 6.7 2.00487 6.70739 1.99934 6.70089 
2.5 6.1 2.47543 6.07171 2.49875 6.10149 
Mean 1.44 7.58 1.44 7.57961 1.44 7.57980 
Variances 0.6830 1.532 0.68219 1.51378 0.68138 1.52834 
Correlation - -0.99617 - -0.99976 - -0.99618 
Third iteration Fourth iteration Fifth iteration 
X3 Y3 X4 Y4 X5 Y5 
0.47039 9.04035 0.50222 9.09583 0.47153 9.03835 
0.81209 8.51668 0.80151 8.49748 0.81283 8.51540 
1.43912 7.55859 1.40009 7.50000 1.43912 7.55851 
2.00420 6.70827 1.99868 6.70179 2.00353 6.70916 
2.47420 6.07322 2.49749 6.10298 2.47298 6.07474 
Mean 1.44 7.57942 1.44 7.57962 1.44 7.57923 
Variances 0.68058 1.51018 0.67977 1.52469 0.67897 1.50657 
Correlation - -0.99976 - -0.99618 - -0.99976 


We have generated the cloned data sets for following nonlinear 
; — 1 _ B pi. ena 
regression models Y = ee Y=A+ T Y = A+ BVX, Y = AX? + 


BX and Y = A + BX + CX? by using the procedure given by [3] and 
presented, respectively in Table 5-9. 
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Table 5. Cloned Data Sets Having the Same Non Linear Regression 


Fit Y = —_: A = 82.97359, B = —36.58871 
A+BX 


Raw data First iteration Second iteration 

X Y Xi Yı X2 Y2 
2.000 0.0615 1.225281 0.102081 1.313470 0.026218 
2.000 0.0527 1.188237 0.102081 1.313470 0.025318 
0.667 0.0334 1.038642 0.017074 0.648053 0.022237 
0.667 0.0258 0.918315 0.017074 0.648053 0.020254 
0.400 0.0138 0.458482 0.014633 0.514769 0.015106 
0.400 0.0258 0.918315 0.014633 0.514769 0.020254 
0.286 0.0129 0.389507 0.013791 0.457862 0.014551 
0.286 0.0183 0.701590 0.013791 0.457862 0.017451 
0.222 0.0083 -0.196641 0.013360 0.425914 0.011090 
0.222 0.0169 0.639830 0.013360 0.425914 0.016789 
0.200 0.0129 0.389507 0.013218 0.414932 0.014551 
0.200 0.0087 -0.121065 0.013218 0.414932 0.011441 

Third iteration Fourth iteration Fifth iteration 

X3 Y3 X4 Y4 Xs Ys 
0.926740 0.028641 0.970763 0.020381 0.777712 0.021073 
0.908248 0.028641 0.970763 0.020104 0.768481 0.021073 
0.833572 0.016874 0.638594 0.019057 0.731203 0.016776 
0.773506 0.016874 0.638594 0.018291 0.701219 0.016776 
0.543963 0.015591 0.572061 0.015855 0.586634 0.016118 
0.773506 0.015591 0.572061 0.018291 0.701219 0.016118 
0.509531 0.015101 0.543653 0.015545 0.569446 0.015852 
0.665320 0.015101 0.543653 0.017056 0.647214 0.015852 
0.216933 0.014839 0.527705 0.013327 0.423385 0.015707 
0.634490 0.014839 0.527705 0.016734 0.631824 0.015707 
0.509531 0.014751 0.522223 0.015545 0.569446 0.015658 
0.254660 0.014751 0.522223 0.013577 0.442217 0.015658 

Table 6. Cloned Data Sets Having the Same Non Linear Regression Fit Y = 


A+- ;A = 0.086506, B = —0.091625 
1+X 


Raw data First iteration Second iteration 
X Y Xı Yı X2 Y2 
2.000 0.0615 2.332049 0.055964 1.805105 0.059008 
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Raw data 


First iteration 


Second iteration 


X 


Y 


Xi 


Yi 


X2 


Yo 


2.000 


0.0527 


1.565860 


0.055964 


1.805105 


0.050796 


0.667 


0.0334 


0.705670 


0.031542 


0.652333 


0.032788 


0.667 


0.0258 


0.506758 


0.031542 


0.652333 


0.025696 


0.400 


0.0138 


0.272456 


0.021059 


0.404582 


0.014499 


0.400 


0.0258 


0.506758 


0.021059 


0.404582 


0.025696 


0.286 


0.0129 


0.257787 


0.015258 


0.296953 


0.013659 


0.286 


0.0183 


0.351251 


0.015258 


0.296953 


0.018698 


0.222 


0.0083 


0.187800 


0.011526 


0.236035 


0.009367 


0.222 


0.0169 


0.325711 


0.011526 


0.236035 


0.017392 


0.200 


0.0129 


0.257787 


0.010151 


0.215011 


0.013659 


0.200 


0.0087 


0.193575 


0.010151 


0.215011 


0.009740 


Third iteration 


Fourth iteration 


Fifth iteration 


X3 


Y3 


X4 


Y4 


X5 


Y5 


2.072216 


0.053842 


1.644782 


0.056682 


1.863838 


0.051862 


1.444278 


0.053842 


1.644782 


0.049020 


1.340783 


0.051862 


0.687722 


0.031054 


0.638878 


0.032217 


0.671312 


0.030599 


0.504364 


0.031054 


0.638878 


0.025600 


0.502137 


0.030599 


0.284091 


0.021273 


0.408885 


0.015152 


0.295141 


0.021472 


0.504364 


0.021273 


0.408885 


0.025600 


0.502137 


0.021472 


0.270142 


0.015859 


0.307342 


0.014368 


0.281892 


0.016421 


0.358695 


0.015859 


0.307342 


0.019070 


0.365715 


0.016421 


0.203334 


0.012377 


0.249424 


0.010363 


0.218200 


0.013172 


0.334572 


0.012377 


0.249424 


0.017851 


0.342948 


0.013172 


0.270142 


0.011095 


0.229361 


0.014368 


0.281892 


0.011975 


0.208864 


0.011095 


0.229361 


0.010711 


0.223486 


0.011975 


Table 7. Cloned Data Sets Having the Same Non Linear Regression Fit Y = 
A + BVX; A = —2.341445, B = 3.011197 


Raw data First iteration Second iteration 
X Y Xi Yı X2 Y2 
1.0 1.1 1.389823 0.669752 1.099958 1.208477 
1.5 1.3 1.536094 1.346503 1.571153 1.390611 
2.0 1.6 1.769221 1.917031 2.033473 1.663810 
2.9 2.0 2.105666 2.419675 2.490122 2.028076 
3.0 2.7 2.764870 2.874101 2.942743 2.665542 

School of Science $ U MT— 13 


Volume 7 Issue 3, 


2023 


Bivariate and Multivariate Data Cloning... 


Raw data First iteration Second iteration 
X Y Xi Yı X2 Y2 
3.5 3.4 3.513706 3.291989 3.392311 3.303008 
4.0 4.1 4.352175 3.680949 3.839463 3.940474 
Third iteration Fourth iteration Fifth iteration 
X3 Ve X4 Y4 X5 Ys 


1.468250 0.816665 1.195129 1.307264 1.541545 0.950454 
1.604771 1.432959 1.637384 1.473126 1.668620 1.511692 
1.820930 1.952519 2.064198 1.721920 1.868668 1.984837 
2.130380 2.410261 2.481144 2.053645 2.153013 2.401687 
2.730323 2.824090 2.891081 2.634163 2.699051 2.778547 
3.404599 3.204646 3.295707 3.214682 3.306735 3.125107 
4.153208 3.558859 3.696127 3.795200 3.976065 3.447676 
Table 8. Cloned Data Sets Having the Same Non Linear Regression Fit Y = 
AX? + BX; A = 0.4,B = 5.0 


Raw data First iteration Second iteration 

X Y Xı Yı X2 Y2 

0 1 0.1933 0.0000 0.0000 1.0107 

1 5 0.9156 5.4708 0.9989 5.0388 
2 12 2.0321 11.7428 1.9970 12.0451 
3 20 3.1495 18.8161 2.9944 20.0064 
4 25 3.7859 26.6907 3.9914 24.9641 
5 36 5.0621 35.3665 4.9881 35.8353 
Third iteration Fourth iteration Fifth iteration 
X3 Y3 X4 Y4 Xs Y5 


0.1932 0.0000 0.0000 1.0211 0.1931 0.0000 

0.9148 5.5123 0.9980 5.0763 0.9140 5.5523 

2.0291 11.7893 1.9941 12.0884 2.0263 11.8338 

3.1436 18.8311 2.9889 20.0114 3.1377 18.8445 

3.7780 26.6377 3.9829 24.9279 3.7701 26.5848 

5.0499 35.2093 4.9762 35.6734 5.0376 35.0546 
Table 9. Cloned Data Sets Having the Same Non Linear Regression Fit Y = 
A + BX + CX?; A = 1.0,B = —0.20,C = 0.20 


Raw data First iteration Second iteration 
X Y Xı Yi X2 Y2 
1.0 1.1 0.99760 1.08571 0.99972 1.09900 
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1.5 1.3 1.52005 1.29286 1.50226 1.30180 
2.0 1.6 1.98929 1.62143 2.00206 1.60299 
2.5 2.0 2.44462 2.07143 2.50109 2.00294 
3.0 2.7 3.05159 2.64286 2.99979 2.70103 
3.5 3.4 3.53894 3.33571 3.49833 3.39798 
4.0 4.1 3.95791 4.15000 3.99676 4.09427 
Third iteration Fourth iteration Fifth iteration 

X3 Y3 X4 Y4 X5 XS 


0.99730 1.08491 0.99945 1.09800 0.99700 1.08412 
1.52239 1.29459 1.50450 1.30362 1.52472 1.29631 
1.99134 1.62433 2.00410 1.60599 1.99338 1.62719 
2.44578 2.07414 2.50217 2.00589 2.44694 2.07680 
3.05120 2.64404 2.99959 2.70207 3.05081 2.64518 
3.53715 3.33401 3.49667 3.39597 3.53537 3.33230 
3.95486 4.14408 3.99355 4.08856 3.95181 4.13819 


4. CLONING FOR MULTIPLE NONLINEAR REGRESSION 
MODEL Y = A(X”) (x) VIA PIVOTS 


In accordance with [4], the current approach was extended to a structure 
of an arbitrary error covariance after discussing data that were independent 
and identically distributed (iid). Let us give the multiple NLRM in (Eq. 4.1). 


Y, = A (XPP) 42, (4.1) 


Where Y is response vector, X = (X 0, X 2) is the covariate data matrix 
(CDM), a = (Ao, 41, A42)" is parameters vector, and E random error 
vector. 

Eq. 4.1 is linearizable due to logarithmic transfomation, then Eq. 4.1 
becomes 1n(Y;) = In(Ay) + Ay In(x{?) + A, m(x®) +E. Setting Y = 
In Y), © = n (xX), ÃO =n (xX), F = In (Ei), Bo =In (40), 
B, = A, and B, = Az we get 
Ÿ = Bo + B® +B, XO +É (4.2) 
where is response vector, = (X: ®) is CDM, B = (Bo, By, Bz)" is 
unknown parameters vector and Ẹ errors vector. When matrix X of rank full 
as column, estimates of £ by ordinary least square (OLS) is b = 
(X°X)-1X°Y, and fitted multiple LRE is 
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Y, = bo + b,x + b,x”) (4.3) 
These consequences follow whether or not Ÿ, X® „and X® are mean 


corrected (MCtD), as here. Due to the MCtD, (Eq. 4.3) , which can be 
written as: 


ĵ = b,x, + boX2 (4.4) 


where P=Y—Y, x =X —XM and x, = O —X@ (In order to 
avoid any loss of generality, the column of 1 in design matrix X is 
eliminated following the imply MCtD because it transforms into a column 
of zeros). 


The identified problem here was to create a new response variable 
vector, Yeione, and a new covariate data matrix, X¢jgne. This can be easily 
accomplished by transposing back to Yejgneand Xcione, Such that 


— fyvet Vv -1t 7 
b= (XcioneXcione) X clone clone 


Alternatively, multivariate CDSs to be required (Yeigne,X es X om e) 
which produced the same multiple NLRM equation as the original dataset 


a) y2) 
(Y, xX®, X2). 


Returning to the case of iid, how generation of CDSs can be 
accomplished via manipulating any one covariate, was exhibited, say xj, 
where j = 1,2., using the steps below. 


1) Iniatially, fit multiple linearizable RM (Eq. 4.4), using MCtD data. 
2) Select a covariate x2. 


3) Let ĵ = k + bax) , where k = b,x, =Y—b2x2. To obtain the 
estimated values of Pg and Xz, perform simple RM of Yg = y — k on x, 
and inverse simple regression of x, on Yg. 


4) Regress x, on x, and obtain X, . Also, obtain x12 = x, — 
Xq(xbx2) txix] = (I — x (xfx) txÉ)xı , where I is the identity 
matrix. 

5) Form Yk, clone 7 Vet > by X12. 

6) On all newly acquired x, and £, , perform multiple LREs of yx cione 
at the same time, where 2, = y(y‘y)~+y‘x, in which y = (1:9) is 
nx 2. 
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7) If you prefer, you can add back Y, ŽO, and X® to the cloned data or 
multiply all of the cloned data by the same scale factor. 


8) Transform back to the cloned data in step 7. 


9) Repeat from one to eight steps to create a series of datasets with identical 
NLRM coefficients. At each iteration, a different possible value of j can 


be chosen. 


Example 4. With uncorrelated data and design matrix of full-rank, consider 
variables X1, X2 and Y in Table 9 and resulting in the multiple NLRM fit 


Y = 1.663079X26163121 y9.2931787 


(4.5) 


For CDSs in Table 10, 1-9 steps specified above were used (first X2 was 
used for manipulation) for which the fitted multiple NLRM equation was 
exactly the same as in (Eq. 4.5). 
Table 10. Cloned Data Sets Having the Same Multiple Nonlinear 


Regression Fit Y = Ao X4 X3? 


Xı X2 Y XA clone X2,clone Y clone 
23.81 11.33 22.76 25.583 10.257 24.988 
15.83 25.92 716.13 37.324 35.179 40.197 

9.46 7.03 8.62 15.945 3.818 16.234 
5.71 29.68 10.98 2.413 14.871 7.852 
85.78 21.81 86.77 49.689 39.486 45.583 
0.37 0.57 0.97 6.672 2.411 4.543 
8.82 11.25 11.82 9.540 9.225 13.577 
8.99 19.01 16.63 5.928 20.744 11.810 
37.65 75.25 67.40 6.780 73.887 19.203 
8.43 8.40 8.81 12.012 4.796 14.364 
16.10 30.30 21.54 6.839 16.206 14.786 
0.64 1.20 1.34 5.718 2.327 5.138 
5.28 6.93 12.38 9.021 22.489 11.379 
30.40 70.18 58.37 5.847 71.736 17.173 
33.66 21.06 29.90 20.152 11.872 25.870 
15.72 11.86 14.54 16.178 6.394 19.093 
8.44 14.53 17.54 7.172 26.020 12.274 
30.20 34.20 29.43 11.444 13.458 21.041 
8.89 8.68 11.41 12.282 8.357 14.702 
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Xi X2 Y X1,clone X2,clone Y clone 
5.14 2.84 5.45 20.372 3.117 14.474 
Figure 2, represent a matrix plot of raw and cloning data in Table 10, 


which show the effect on X2 done by orthogonal manipulation as described 
in steps 1-9 of the algorithm. Bivariate relationship strength between X2,clone 
and Ycione is much weaker than X2 and Y. However, this is not the case with 
Xi clone and Yclone , because the manipulation was not done with X1. 


Cloned Data for Multiple Nonlinear Regression 


20 40 60 80 


0 


20 40 60 80 


0 


Matrix plot of raw data (X,, Xz, Y) and cloned data (Xj iene. X2cione: Y cione) 


Figure 2. Matrix Plot of raw and cloned data 
5. DISCUSSION 


This study showed that the parameter estimates of the original datasets 
discussed in this artical and their generated cloned datasets were identical. 
As a result, it was identified that data cloning had the potential to be used 
in a wide range of applications, including data encryption, visualization, and 
smoothing. The application of encryption was particularly intriguing 
because it can be used to generalize the databases even when regression 
modeling was not desired. In prior literature, cloned datasets were 
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generated for linear regression models. However, it had equal importance 
to be generated for the nonlinear regression models. In this context, new 
methods can be developed for nonlinear regression models to conduct 
cloning for the datasets or databases. 


5.1. Conclusion 


CDSs have been presented for bivariate and multivariate NLRMs that 
have linearizable regression functions including 


AX®, ABX, Ae®X,—"_, A+ Z,A + BVX, AX? + BX, A + BX + CX? 


and Ag)" Go with exactly the same nonlinear regression 
coefficients. In terms of bivariate LRFs, the response and a covariate of the 
CDSs collapsed to their means, which had smaller variability when 
compared to the original dataset. 
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ABSTRACT 


Water is the basic need for the survival of all living organisms; however, its 
contamination has caused disastrous effects on human health. Almost 80% 
of the total population in Pakistan is forced to use unsafe drinking water due 
to the scarcity of clean water sources. To assess, the risk involved and to 
monitor the amount of Arsenic and bacterial growth in drinking water, water 
samples were collected from different government schools and 
underdeveloped areas of Multan, Pakistan. Moreover, water samples were 
collected by randomly visiting areas and were sent for laboratory analysis 
to determine and detect the arsenic amount along with other parameters 
including aesthetic, physical, chemical, and bacteriological growth. The 
results indicated that all the data collected from Govt. schools and 
developed locality were having a large amount of arsenic in drinking water. 
A huge amount of Coliform growth was found in the Peer colony water, 
which indicated a highly unhealthy intensity of pollutant water for the 
population residing in this area. Therefore, urgent and regular monitoring 
of drinking water in Multan was required to prevent such contamination. 


Keywords: arsenic, coliform, health risk, water pollution, water scarcity 
1. INTRODUCTION 


Water holds the second most significant position after oxygen, which is 
essentially important for human survival. Therefore, water is an essential 
component of life not only for human beings but for all living organisms. 
In some living organisms, approximately 90% of the body weight consists 
of water. The majority of the world's population is suffering from water 
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scarcity, due to unsafe and unsuitable water consumption, which has been 
regarded as a fundamental right for nations. As indicated in a study, the 
heart and brain mainly comprise almost 73% of water, whereas 83% of 
water is found in the lungs. Water plays a vital role in our body by 
maintaining the body temperature, functions as a shock absorber for 
sensitive parts of the body like the brain and spinal cord, which aids in the 
growth of body cells, forms saliva that helps in digestion, and maintains the 
moisture in the mucosal membrane. Various human activities have caused 
impurity and toxicity in water, which has resulted in water contamination, 
which is a leading health problem nowadays. Harmful metals, such as 
arsenic, iron, and cobalt are easily and excessively found in water. These 
metals have affected different organs of the body, such as kidney 
dysfunction, cancer, nervous system disorder, skin lesions, and immune and 
gastrointestinal disorders. According to a case study, almost five million 
children have lost their life because they were subjected to bad water 
consumption and toxic water intake [1]. 


In many areas, fecal matter is mixed with water due to poor sanitation 
systems and leakage of pipes in poorly developed areas. This contamination 
may worsen health problems, which may lead to serious health problems 
[2]. Different pathogenic bacteria are found in drinking water such as 
Shigella sonnei, Salmonella enterica, and other bacteria. These bacteria 
were chiefly disseminated by fecal excretion [3]. Coliform bacterium was 
found in the feces of warm-blooded animals, which cause severe disease 
when exposed to water [4]. 


However, Pakistan has faced water toxicity issues due to its poor 
sanitation systems. In some areas, these detrimental chemicals and 
pathogens have caused major health problems by putting massive people’s 
life at stake. According to WHO and UNICEF, approximately 1.1 billion 
people worldwide are suffering due to the nonexistence of pure drinking 
water [5]. Many other water-borne diseases like hepatitis, diarrhea, and 
cryptosporidium are increasing rapidly due to the flawed sanitary system. 
As per Child Rights International Network, nearly 250,000 children have 
faced a high death rate, due to the high impurity and toxicity of drinking 
water [6]. In industrial areas of Pakistan, water gets contaminated when 
exposed to injurious chemicals. Furthermore, it was also identified that in 
underdeveloped areas of Pakistan, when drinking water is mixed with 
sewage water it caused many health issues even leading to death [7]. The 
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quantity of arsenic found in drinking water was quite high in secondary 
schools of Vehari (mostly ranging between 12.9ug/L-19.5ug/L), which 
predominantly affected the health of students [8]. Highly populated cities 
like Karachi, Lahore, Rawalpindi, and others face severe drinking water 
pollution due to human activities, rendering it unsuitable for consumption 


[9]. 


Many developing countries like Pakistan are facing the common issue 
of water scarcity. Pakistan has ranked as the 7" country, which faces water 
pollution and water scarcity problem. Exterior as well as interior water is 
contaminated with arsenic. Arsenic is a carcinogen that causes cancer [10]. 
Chlorination is an effective method to extract toxic chemicals and microbes 
from the water, making it suitable for drinking [11]. 


2. MATERIALS AND METHODS 


The current study was conducted in different government schools and 
some local underdeveloped areas of Multan like Mumtazabad, Peer Colony, 
Chungi no. 8, Northern Bypass, and Ghanta Ghar Multan, to find out the 
amount of arsenic in the drinking water. The Environmental Protection 
Agency (EPA) has set the standard for arsenic in drinking water to be 10 
micrograms per liter. This study considered all the physical, bacteriological, 
synthetic, and aesthetic analysis reports issued by the water quality 
laboratory in Multan to conduct further analysis. 


2.1. Sample Collection 


For this purpose, the researcher randomly visited government schools 
and collected water samples from the cooler, electric cooler, drinking water 
bottles, water fountains, and taps. Sterile 250mL plastic bottles were used 
for the sampling purpose. During collection of samples gloves were used by 
the collector. The bottles were filled up to 240mL leaving the headspace of 
the bottle, capped the bottles immediately, and delivered to the lab on the 
same day of their collection. 


2.2. Bacterial Detection 


For the detection of total coliform, fecal coliform, and E. coli petri, 
petrifilm calculation plates were used. Place Petrifilm EL plate on a leveled 
surface and lift the top to place 1ml of sample onto the center of the bottom 
film. Now release the top film and allow it to drop. Place the spreader on 
top of the film and gently apply pressure to distribute it evenly. Lift the 
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spreader and wait for one minute for the gel to solidify. Now incubate the 
petri film count plates on the clear side up at a temperature of 30°C for 
almost 48 hours. In, the end the bacterial colonies growing on petrifilm were 
identified with a magnifier. 


2.3. Arsenic Test 
Arsenic was analyzed in the drinking water samples by using a hydride 


generation atomic absorption spectrometer (HG-AAS). All samples were 
tested thrice for better results. 


3. RESULTS AND DISCUSSION 


Drinking water samples from different schools identified a high level of 
arsenic (15ug/l-75ug/l). None of the selected schools’ arsenic level in 
drinking water was normal. The graph given below presented a level of 
arsenic in school water available for school children. 


Arcenic level 


S1 S2 s3 S4 s5 S6 S7 s8 s9 S10 


=m Arsenic (ug/l) recommended limit of arsenic (ug/l) 


Figure 1. Arsenic Level in Drinking Water of Different Schools 


The first step was to collect samples from different schools located in 
different areas of Multan. The above graph presents the level of arsenic in 
different schools. The schools are presented with symbols and numbers 
from S1-S10. According to the results sample collected from S5 indicated 
arsenic 15ug/l, which is the lowest value as compared to other samples 
collected, however, they were still unsafe to drink. The highest value of 
arsenic was evaluated in S8 School with 70ug/l. These values of arsenic 
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were also informed in schools and later schools set up filter water for 
children. 


A direct link between drinking water and health was noticed in the 
samples. The consumption of arsenic in drinking water can lead to various 
problems like skin diseases, stomach diseases, and cancers. Also, in order 
to check the concentration of arsenic in water, its distribution was also 
analyzed across different locations in Multan City. 


Arcenic level in Multan 


Mumtazabad peer colony Northern bypass Chungi # 08 Ghantaghar 


=m Arsenic (ug/l) ——recommended limit of arsenic (ug/l) 


Figure 2. Arsenic Data in Drinking Water of Different Locations in 
Multan 


The above graph presents arsenic values tested from the water samples 
collected from different areas of the Multan, Pakistan. These areas were 
selected randomly. The lowest level of arsenic found in selected areas was 
Mumtazabad with 25ug/l and highest value in areas Northern Bypass and 
Chungi no.8 with 75ug/l. However, even the lowest value found in 
Mumtazabad water was still high to consume. The other two areas Ghanta 
Ghar and Peer colony both were found with 50ug/l arsenic in drinking 
water. The presence of arsenic in these areas can be due to high amount 
present in the ground water plus the rusting of pipes and the contamination 
of water due to high industrialization in these areas. A small survey was 
also conducted related to different skin, eyes, kidney, stomach, and other 
diseases in order to recognize conditions of people using this water. This 
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area according to our small survey indicated people with skin, stomach 
diseases, and eyes issues each house. Although detail planning and data 
collection is required to find the complete detail of types of skin, stomach, 
and eyes diseases related to arsenic in water. All the data collected was 
depicted in the graphs to determine the purity of water. In drinking water 
supplies, arsenic is a significant problem because of its toxicity. In this 
research, the aesthetic, physical, chemical, and bacteriological analysis 
reports were incorporated, which were issued by the water laboratory, 
Multan, Pakistan. 


According to previous research, arsenic was calculated (0->50ug/1) in 
2005 Multan City, while [12] the updated data of 2022 presented the 
maximum value calculated of arsenic in the Multan, which is 75ug/1. There 
are limitations applied here but still it gave us an idea about the increase of 
arsenic level in the city water from 2005-2022, which is an alarming 
situation. 


There is another issue of bacterial colonies in the drinking water. 
Various samples from different areas of Multan for E. coli and Coliform 
bacterial colonies were tested. None of the water was found with disease 
causing bacteria except Peer colony water. The petri film was observed very 
carefully and was found with colonies of total coliform. Five colonies were 
found countable, while others were not in range of counting due to the 
presence on gridline and out of the circle. The coliform bacteria present in 
the water sample of Peer colony was further evaluated for the presence of 
fecal coliform. Peer colony was found with 2 colonies of fecal coliform. 
These results were evaluated thrice and calculated as per ml of sample and 
the presented colony numbers are the mean of all results. The results 
presented in Table 1 indicates absence of E. coli and presence of coliform 
bacteria with specifically indication of some colonies of fecal coliform 
bacteria. Such type of bacteria was hazardous for human health. This 
presence opens a new research study to be setup in Peer colony. The 
presence of total coliform and fecal coliform bacteria in water indicated risk 
of contracting a water borne illness increase in this area. There are other 
publications which have presented a high percentage of arsenic in Multan 
but no specific areas were mentioned that can be compared with future 
studies [13]. In parallel, another ongoing study sheds light on arsenic levels, 
underscoring groundwater contamination and echoing our findings of 
elevated concentrations. Together, these revelations emphasize an 
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immediate call for proactive mitigation strategies to address arsenic 
contamination in Sindh's water sources. Elevated levels, influenced by 
natural processes and human activities, exceed WHO limits, necessitating 
urgent action to mitigate health risks and ensure the safety of the water 
supply [14, 15]. 


Table 1. Identification of Bacterial Colonies in the Water of Peer Colony 
Houses (Per Sample) 


Sr Heights Desirable 


No. Water Quality Parameters Leva Result 
1 Total Coliform 0/100 ml (PSQCA) 5 
2 Fecal Coliform 0/100 ml (PSQCA) 2 
3 E. coli Negative Negative 


A prior study conducted on the localized areas of Multan presented no 
signs of coliform bacteria, however, this study found the presence of total 
coliform and fecal coliform in Multan city from every one out of five areas 
of the collected samples. This difference can be due to the area selection or 
it can be due to clean water present in 2005. After the positive result of peer 
colony bacterial contamination, there should be a selection of many more 
areas to confirm other species and chemicals present in the water. 


3.1. Conclusion 


The current study identified that in many cases, groundwater was found to 
be polluted with high amounts of arsenic, which has caused major health 
issues, especially in young children and adults. Notably, if groundwater is 
recharged directly from the land surface, which has a historical association 
with industrial usage or the disposal waste, there exists a significant risk of 
high arsenic concentrations in the groundwater. For this purpose, new 
technologies for purification should be developed to achieve high quality 
drinking water to avoid any future health disasters and exposure due to the 
toxic elements present in the drinking water. Hence, it was identified that 
Multan water still requires an extensive study with large study samples to 
cover major external and internal factors in order to find specific areas 
where arsenic level is high, which should be considered serious for health 
issues. 
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ABSTRACT 


Breast cancer is a significant health problem for women globally; 
however, timely detection can reduce female morbidity and mortality. 
Early breast screening has become imperative for all women, though, 
adequate screening facilities are necessarily required in developing 
countries like Pakistan, where breast cancer is a leading cause of death. To 
encounter this chronic disease, various image processing techniques have 
been introduced to automatically diagnose breast cancer from digital 
mammograms. The current study deployed data from a population of 35 
participants. The mammograms used for screening were 5 normal, 15 
benign, and 15 malignant patients. The breast images were marked by the 
radiologist and the system was trained with normal, benign, and malignant 
classes. Moreover, Multilayer Neural Networks (MNN) based texture 
analysis methodology was adopted to distinguish normal, benign, and 
malignant breast images. Reportedly, an automated approach was used to 
detect breast conditions after conducting the analysis of digital 
mammograms. Statistical parameters, namely sum, mean, variance, 
standard deviation, kurtosis, skewness, energy, and entropy were 
calculated, analyzed, and compared for the normal, malignant, and benign 
breast images. The results indicated a 100% accuracy after the analysis. 
The results of the extracted statistical parameters were promising and 
reliable in distinguishing between normal, malignant, and benign breast 
mammograms, again indicating the need for early detection of the disease 
to minimize the risk of breast cancer among women. 


Keywords: breast cancer, benign and malignant cancer, mammography, 
MNN, neural network 
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1. INTRODUCTION 


Breast cancer is one of the most common cancers among women, 
which affects women's physical and mental health and is even life- 
threatening globally. The overall frequency is increasing, and the risk 
factors are well-documented [1]. 


There are several features, which are responsible for this disease like 
obesity, radiation exposure, family history, fat-enriched food, and pre- 
malignant injuries. Breast cancer once detected can become worse with a 
change in lifestyle [2, 3]. Breast cancer is the second major cause of 
mortality worldwide and in Pakistan every 1 woman out of 9 is suffering 
from this chronic disease. This can be due to the lack of medical facilities 
and unawareness, which are major causes of delayed diagnosis [4, 5]. It 
has been a challenging task to detect breast cancer at an earlier stage to 
lower the mortality rate. The risk of breast cancer among women aged less 
than 40 years is smaller, whereas an increased rate has been noticed in 
women between 45-55 years, resultantly causing an increase in mortality 
rate [6]. Mammography is an affordable, low-radiation test recommended 
for early detection of breast cancer [7, 8]. This test can be performed on 
both pre-symptomatic and symptomatic women; however, it is less 
efficient for women who are lower than 40 due to dense breasts [9]. 
Though, women at later stages can feel a lump, a hard loop within the 
breast or under the arm, swelling, changes in the mass of the breast, 
redness, and dimpling or creasing of skin. Countries with low-income and 
underdevelopment, such as Pakistan, face a scarcity of resources for 
essential diagnosis or detection. Delays in medical care may contribute to 
the spread of breast cancer’s morbidity and mortality [10]. Michell et al. 
Identified that the addition of digital breast tomosynthesis (DBT) in 
conventional mammography increases the accuracy rate in cancer 
assessment [11]. 


The human nervous system is a compound neural system and our body 
collects external signals through dendrites. These signals are collected till 
a threshold level to make an action. A large number of neurons are 
interrelated to study input and optimize its output by means of computing 
and processing. Artificial Neural Network (ANN) is used in training a 
learned system by processing these signals in medical imaging so, that 
decisions can be automated. Computer-aided Detection Systems (CAD) 
trained using machine learning models, help radiologists, medical 
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practitioners, radiation oncologists or physicians in the accurate and more 
definite detection of breast cancer in the early stage [12]. CAD systems 
have been used in many medical institutes for clinical diagnosis as 
referred by doctors [13]. For the feature extraction process, curvelet 
moments are also very beneficial in giving reduced features and enhancing 
accuracy in the diagnosis of breast cancer [14]. An automated feature 
extraction process is utilized in diagnosing architecture distortion in the 
breast, which causes malignancy by applying texture models and vector 
machine classifiers on the region of interest [15]. A multivariate feature 
can also be used for extraction by applying a genetic algorithm and 
obtaining feature descriptive in two categories, which classify malignant 
and benign breasts and reduces the workload of radiologists [16]. 
Furthermore, the grey level co-occurrence matrix is used for feature 
extraction in different orientations to differentiate masses and non-masses, 
by combing the grey level with the wavelet approach through which the 
results are enhanced [17]. ANN-PSO (artificial neural network with 
particle swarm optimization) method has been used by Zhang et al to 
extract the intensity and shape features, which tested and then 
distinguished normal and abnormal breast mammograms [18]. The particle 
swarm optimized wavelet neural network (PSOWNN) technique was 
applied to breast mammograms by Sathyaa et al and the Receiver 
Operating Characteristic (ROC) curve indicated the sensitivity and 
specificity as indicated by Dheeba et al. [19, 20]. 


From the analysis of texture, an automated method was proposed to 
classify different breast tissues as normal, malignant, or benign. This 
method provided the radiographic details of breast lesions. 


Breast tumors can be benign or malignant. Benign is not dangerous, it 
can grow but doesn’t spread in surroundings. However, malignant is 
cancerous and can spread in the entire body and cause destruction [21, 22]. 
Therefore, in this study, an automated approach has been proposed to 
classify breast cancer as normal, benign, or malignant. This work explored 
ways to prevent this disease and provide a new taxonomy for breast cancer 
risk assessment. 


2. MATERIALS AND METHODS 


A dataset of mammograms, correspondingly collected from 35 
subjects was obtained from the Bhawal Victoria Hospital. The women 
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who were considered for scanning ranged from age 45-55 years having 
either normal, benign, or malignant breast cancer stage. All patients were 
duly informed and the ethics review board signed for further process. 
Medial-lateral oblique and craniocaudal images were acquired and 
statistical mammographic data were analyzed. The following steps were 
taken to successfully classify breast conditions including normal, benign 
or malignant. 


i. The mammograms having high intensity and contrasts were input to 
the developed system. The developed system consisted of an 
architecture based upon 12-layered multi-layered neural networks 
(MNN). 


ii. Segmentation of Region of Interest (ROI), which had a suspicious area 
with distinguished surroundings was marked by the domain experts. 


iii. Texture features, namely, sum, mean, variance, standard deviation, 
kurtosis, skewness, energy, and entropy extracted from each of the 
segmented ROIs belonging to normal, benign, and malignant stages. 


iv. Training of the system using these extracted texture features through 
MNN using Matlab2020a was conducted on 70% of the total data. 


v. Evaluation of the proposed developed system was based on the 
standard evaluation measurement, such as accuracy. 


3. RESULTS 


In this research, statistical methods were used to extract texture 
features. Using image classification and segmentation methods was a more 
reliable approach, which ensured higher reliability in the obtained results. 
The feature extraction method underlying the statistical procedure 
consisted of the sum of the following parameters, namely mean (pixel 
value), variance (square root of deviation), energy (uniformity), standard 
deviation (pixel deviation from the mean), skewness (pixel asymmetry), 
kurtosis (flatness or peak), and entropy (random measurements). Normal, 
malignant, and benign breast mammograms and their bilateral statistical 
data were analyzed. There are usually two models for image projection, 
namely Craniocaudal (CC), a view from the apex, which allows for 
enhanced imaging of the central and internal breast regions, and 
Mediolateral Oblique (MLO) that provides an improved view of the gland 
from a specific side view. 
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(b) (d) 
Figure 1. (a) Case 1 (L-MLO), (b) Case 2 (R-CC), (c) Case 35 (R-MLO), 
and (d) Case 35 (R-CC) 


(a) (b) (c) (d) 
Figure 2. (a) Case 7 (R-MLO), (b) Case 7 (L-CC), (c) Case 19 (L-MLO), 
and (d) Case 19 (R-CC) 


Figures 1 and 2 show a normal case in which no tumor was detected 
on mammography screening of the two selected subjects. 


(a) (b) (c) @) 
Figure 3. (a) Case 1 (L-MLO), (b) Case 1 (L-CC), (c) Case 2 (L-MLO), 
and (d) Case 2 (L-CC) 


Figure 3 contains randomly selected left mammograms from two 
different subjects. The regions highlighted in these MLO and CC images 
are malignant breast lesions. Benign areas are brighter due to their higher 
intensity. Figure 4, shows right mammograms of two cases containing 
malignant areas in MLO and CC image projections. In Figure 4, Case 3 R 
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(MLO and CC) shows abnormalities, lesion-like masses, foci, asymmetric 
dense margins, and unclear partially malignant pathology (a, b). Case 4 R 
(MLO and CC) showed abnormalities, with lesion-like mass morphology 
and suspiciously malignant irregular margins in parts (c, d). The selected 
area is the tumor-containing lesion. 


(a) (b) (c) (d) 
Figure 4 (a) Case 3 (R-MLO), (b) Case 3 (R-CC), (c) Case 4 (R-MLO), 
and (d) Case 4 (R-CC) 


A 2 3 


(a) O) (©) (d) 
Figure 5. (a) Case 1 (L-MLO), (b) Case 1 (L-CC), (c) Case 2 (L-MLO), 
and (d) Case 2 (L-CC) 


(a) (b) (c) (d) 
Figure 6. (a) Case 3 (R-MLO), (b) Case 3 (R-CC), (c) Case 4 (R-MLO), 
and (d) Case 4 (R-CC). 

Mammograms of the left breast for case 1 and case 2 are shown in 
Figure 5, with the mediolateral and craniocaudally view with affected 
highlighted areas. This study has dealt with images of benign breasts 
having calcification, pleomorphic distribution, and architectural distortions 
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as shown in Figure 6. Case 3 R (MLO) is abnormal. Lesion-type 
classification, pleomorphic distribution, and clustered benign 
mammogram are shown in Figure 6 (a). Case 3 R (CC) is abnormal lesion- 
type calcification, pleomorphic distribution clustered, and architectural 
distortion margins speculated benign as shown in Figure 6 (b). Case 4 R 
(MLO and CC) shows abnormality lesion-type mass and shape round 
margins obscured benign as shown in Figure 6 (c, d). Red enclosed areas 
in the normal and benign breasts are the regions of interest, which are 
segmented and processed for statistical feature extraction of 
mammograms. 


Texture parameters extracted from normal ROIs of the selected 35 
patients age ranging from 45-55 are shown in Table I. Texture analysis 
was performed on 35 women with bilaterally normal breasts. Data from 
all the patients were collected and summarized below, showing medial- 
lateral oblique and craniocaudal views. 


Table 1. Statistical Parameters in Mediolateral Oblique and Craniocaudal 
view of Normal Right and Left Breast 


Normal (MLO) Normal (CC) 
aaa (Average values) (Average values) 
Right Breast Left Breast Right Breast Left Breast 
Sum 53.0 50.00 43.0 42.00 
Mean 81.3 81.67647 66.4 65.470588 
‘Variance 68.0 6768.161 63.0 6470.7764 
~ 82.0 82.04567 79.0 80.050905 
Skewness 0.49 0.455200 0.86 0.8279944 
Kurtosis -1.20 -1.272985 -0.60 -0.7750828 
‘Entropy 26.0 2651414 21.0 211126411 
Flatness 0.46 0.475726 0.52 0.53155200 
Energy 3.2 E-05 3.268E-05 4.0 E-05 4.16294E-05 
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Table 1 provides textural parameters, such as, sum, mean, variance, 
standard deviation, skewness, kurtosis, entropy, flatness, and energy, 
along with maximum, minimum, and mean ranges for all normal 
mammograms. 


The statistical characteristics of all the patients with malignancies were 
determined and presented in Table 2. 


Table 2. Numerical Analysis of Right and Left Malignant Breast in MLO 
and CC view 


oa Malignant (MLO) Malignant (CC) 
Statistical (Average values) (Average values) 
Parameters : s 
Right Breast Left Breast Right Breast Left Breast 
Sum 5754958.9 5665979.1 4779375.86 4332792.4 
Mean 78.371429 77.400000 62.6571429 59.828571 
Variance 6171.9399 6381.0684 5749.74924 5489.5443 
pranan 78.195735 79.510718 75.2393183 73.539115 
deviation 
Skewness 0.5173302 0.5314520 0.89566069 0.9520324 
Kurtosis -1.0368810  -1.0458932 -0.48604212  -0.2584942 
Entropy 30025540 28286342 23421773.6 21287130 
Flatness 0.4867556 0.4924317 0.55017747 0.5496571 
Energy 2.727E-05 2.9818E-05  3.49926E-05  3.8034E-05 


Table 3 summarizes the texture parameters of the image window 
pixels, including the maximum, minimum, and mean values of the left and 
right benign breasts, with promising results. Reportedly, maximum and 
minimum values are now equal to absolute values. 


Table 3. Pixel Parameters of Right and Left Benign Breast Images 
Benign (MLO) 


Benign (CC) 


Statistical (Average values) (Average values) 
Parameters : A 
Right Breast Left Breast Right Breast Left Breast 
Sum 5798909.86 5172456 4832506.3 5066756.51 
Mean 73.2285714 66.485714 62.200000 65.2857143 
Variance 5685.44068 5397.0492 5488.9733 5285.80702 
Standard 74.9817381 73.053862 73.571646  72.37847995 
deviation 
a Scientific Inquiry and Review 
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Statistical Benign (MLO) Benign (CC) 
ee (Average values) (Average values) 
Right Breast Left Breast Right Breast Left Breast 

Skewness 0.68503344  0.7359084 0.9761847 0.80500920 
Kurtosis -0.7213016 -0.604295 -0.1948880 — -0.41404877 
Entropy 28338321.9 25325634 214942223 24647141.59 
Flatness 0.51946173  0.5318824 0.5694440  0.536184773 
Energy 2.727E-05 2.9818E-05 3.49926E-05 3.8034E-05 


The results for the left and right benign breasts are shown in Table 3. 
To obtain reliable results and conclusions about the texture parameters of 
all normal, malignant, and benign breast images, comparisons were made 
numerically, to classify the results of normal and benign breast data by 
using MNN. The proposed approach was analyzed and found to be a 
reliable approach for the assessment. Digital data from previous 
mammograms of benign and healthy breasts were used for differentiation. 
Texture parameters are compared in Table 4 for maximum and minimum 
values for both of the selected cases. 


Table 4. Comparison of Right MLO view Texture Parameters of Normal 
and Benign Breast Mammograms 


be Right normal Right benign 
Diaa {MLO} (MLO) % Differences 
(Average Values) (Average values) 

Sum 5383731.714 5798909.86 -7.1595896 
Mean 81.34285714 73.2285714 11.0807648 
Variance 6807.459804 5685.44068 19.7349544 
onda 82.28868257 74.9817381 9.74496545 
deviation 

Skewness 0.491937579 0.68503344 -28.187801 
Kurtosis -1.20984084 -0.7213016 67.730231 
Entropy 26957357.8 28338321.9 -4.8731329 
Flatness 0.468915193 0.51946173 -9.730559 
Energy 3.23365E-05 2.75458E-05 17.3917621 


4. DISCUSSION 


Normal, benign, and malignant mammograms were analyzed and the 


extracted textural 
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craniocaudal (CC) views. Subtle and obscure areas are enhanced using 
algorithms to make areas of interest easier to see and texture extraction 
methods are more useful in these areas. In this research, mammograms 
were used to train and test the proposed model to automatically classify 
normal, benign, and malignant breast images based on their texture. Each 
mammogram contains 20 slices of both views (MLO and CC) along the 
left and right breast. Therefore, 35 x 2 x 20 x 2 slices were used for 
training and testing the reported system. Of these total 2800 slices, 1848 
slices were used for training the reported system and the remaining 952 
slices were used for testing the reported system. Reportedly, two clues to 
malignant tumors are identified, namely microcalcifications and masses. 
These texture features were used to train the proposed method using 
MNNs and then compared to the literature [23, 24], both types yielded up 
to 90% of accuracy in discriminating benign or malignant abnormalities. 


The statistical analysis examined each category's sum, mean, variance, 
entropy, skewness, and other pixel features of left and right breasts. 
Variations in these parameters, such as maximum, minimum, and mean 
values, were used to classify normal and benign tumor, including breast 
and the areas of microcalcification. In this work, MLO-mapped 
parameters of mammograms of healthy and malignant breasts were used in 
which clear differences in all values were noticed. 


Differences between all parameters and maximum, minimum, and 
mean values demonstrated the intelligent classification of normal and 
benign breasts. This study’s results were promising for various statistical 
approaches to describe and analyze image textures. By analyzing the 
variance of all outcomes, it was concluded that a distinction is required 
between the normal and benign cases with 100% accuracy. Different 
values are also agreed upon for normal and abnormal benign images. The 
observed variance for normal cases ranges from 8831.9-3896.8 and all 
findings for abnormal benign cases are above the maximum for normal 
images. Therefore, the variance parameter gave 100% results from benign 
to normal images. Furthermore, it was observed that normal and abnormal 
benign cases could achieve different ranges of kurtosis 88% of the time, 
with 89% of perfect results. For the rest of the parameters, a clear result 
was also noticed. This demonstrated good results obtained in a large 
number of breast cancer cases and benign cases. These were further 
compared on various parameters of normal, benign, and malignant breasts 
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by plotting kurtosis, energy, flatness and distortion on the x-axis and their 
intensity on the y-axis. Data were obtained from random patients and 
parameters were plotted with intensity. The variation in values of normal, 
malignant, and benign breast mammograms can be clearly observed 
graphically. 


Figure 7. (a) The normal breast parameters of right and left mammogram 
in MLO and CC view, (b) Parameters of benign breast images along with 
intensities, and (c) Malignant breast mammogram parameters along with 
intensities. 


Figure 7(a) shows the changes in kurtosis, energy, flatness, and 
skewness of left and right mammograms of a normal breast in 
craniocaudal and medial-lateral views. The parameters compared with the 
benign mammography parameters are shown in Fig. 7(b). In Fig. 7(b), the 
hyperintensities of the right benign image in CC view and MLO show the 
variation of values in normal parameters. The kurtosis parameter has a 
high intensity, which decreased at the energy point and increased again 
with flatness and strain but increase in the value at cc and then decreases 
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gradually. Comparing these two charts with the malignant breast image, 
the deviations from normal values are larger than the benign image 
parameters, as shown in Fig. 7(c). Intensity fluctuations in right malignant 
CC views are rather instantaneous and different from benign and normal 
breast parameters. The left malignant tumor (MLO, CC) and right 
malignant MLO also showed strange curves, which were different from 
the normal breast curve. This artificial texture analysis approach would 
play an important role in the future. Furthermore, the research aims to test 
these methods with other classifiers or segmentation schemes and in more 
realistic texture interpretation situations. These results showed that there is 
considerable performance variability between different texturing methods. 
These results are supported by previous comparative studies in this field. 
Hence, the feature extraction results may assist radiologists in the early 
screening of breast mammograms and discrimination between healthy and 
unhealthy breasts by analyzing the numerical data. 


4.1. Conclusion 


The current study, aimed to propose a multi-layered neural network 
(MNN) approach to demonstrate the potential to distinguish between 
healthy, benign, and malignant breast mammograms in breast region 
analysis and segmentation. Texture analysis and comparison parameters of 
statistical data such as mammogram sum, variance, mean, kurtosis, 
standard deviation, skewness, energy, entropy, and flatness were used to 
obtain the results. It was obvious from the results that this feature 
extraction modality was reliable for diagnosing cancerous or malignant 
breasts. Therefore, it was recommended that early detection could increase 
the survival rates among women possessing breast cancer. Moreover, 
screening tests are mandatory for all women, especially in developing 
countries like Pakistan, where cancer mortality is relatively higher than in 
other developed countries. 
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ABSTRACT 


A microelectromechanical system (MEMS) is a diminutive machine having 
electronic and mechanical components with a size ranging from 20 um-1 
mm. In this present-day world, MEMS fabrication techniques have remod- 
eled the conventional approaches towards system fabrication. Microfluidics 
is an eminent domain of MEMS in which small volumes of fluids are disci- 
plined in micro-channels having dimensions in the submillimeter to achieve 
the desired outputs. Microfluidics have revolutionized the realm of compact 
system fabrication through preeminent inventions like lab-on-a-chip tech- 
nology. Microchannels of various architectures are fabricated to employ 
microfluidic systems depending upon the required function of the device. 
In ophthalmology, Retinal Vein Occlusion (RVO) is an ailment in which 
small veins that take away blood from the human eye's retina are blocked 
or fissured, causing vision loss. Therefore, in this study, four micro-chan- 
nels with different architectures, namely, sinusoidal, U-shaped, spiral, and 
curvilinear, were simulated by using the fuzzy technique to investigate the 
optimization of fluids for the implantation process to fix the RVO elixir. 
The two most critical parameters in retinal vein flow rate and velocity were 
taken at the output for optimization. Hence, fuzzy fluidic simulation re- 
vealed that curvilinear micro-channels were the best fit for biomedical im- 
plantation to treat RVO malady. 
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1. INTRODUCTION 


A device with electrical and mechanical components as its integral parts 
and a size less than 1mm but greater than one micron falls in the category 
of MEMS [1]. Applications of MEMS have greatly revolutionized almost 
every field of scientific research. Additionally, MEMS based devices are 
widely employed in routine life; including automobile airbag accelerome- 
ters, keyless entry mechanisms, heat exchanging systems for cooling mi- 
crosystems, lab-on-a-chip technology, micro drug delivery mechanisms, 
and many more. In the 1990s, an essential sector of MEMS, microfluidics 
study, came into the limelight with its astonishing applications in the fabri- 
cation of MEMS-based devices. In microfluidics systems, small volumes of 
fluids, usually the norm in the range of microliters (ul) to picoliters (pl) are 
controlled and manipulated in micro-ducts having dimensions from tens to 
hundreds of micrometers to obtain desired outputs [2]. Inkjet printers 
proved to be a prominent commercial success for microfluidics-based 
MEMS [3]. Microfluidics-based devices greatly benefited medical science 
by providing biochips, especially diagnosis applications [4]. 


Micro-channels are an imperative part of any microfluidic device [5]. 
The small volume of fluids is contained and manipulated in micro-channels 
of hydraulic diameter customarily below thanlmm. Originally the concept 
of micro-channels was suggested by Tuckerman and Pease of Stanford 
Electronics Laboratories [6]. Micro-heat exchangers have widely employed 
micro-channels containing small qualities of fluids for heat exchange. These 
micro-heat exchangers play a significant role in the gas turbine engines, air 
conditioning, and heat pumps. It is because of the modest velocities at which 
fluids moving via microchannels, low Reynolds numbers are obtained [7]. 
Primarily, this paper is based on Retinal Vein Occlusion (RVO) for all three 
micro-channels. The human eye's lens focuses on all the incoming light onto 
the retina, located in the rear of the eyeball. The optical nerve transmits 
information from the retina to the brain, where it is processed to recognize 
images. Retinal vein occlusion (RVO) occurs when a clot blocks blood flow 
in the retinal vein. Hence, this research describes the maximal and typical 
blood flow rates in the human retinal veins. 
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1.1. Microchannel Shapes 


Depending upon the intended purpose, micro-channels of different mor- 
phology were concocted. For example, to reinforce the flow pattern like 
parallel, serpentine, crossed, oblique, and ratchet morphology of micro- 
channels could have adopted the architecture application. The shape of the 
microchannel is always an important factor to obtain the required results. 
Computational fluid dynamics (CFD) is often employed to optimize the de- 
sired microchannel shape [8]. The design of different alphabets like J, Y, 
T, U, and S can be a choice for purpose-based on its applications [9]. The 
other viable architectures include triangular, spiral, double spiral, trapezoi- 
dal, rectangular, circular, square, coil, double coil, curved, curvilinear, con- 
centric, and sinusoidal with ascending and descending orders. The preferred 
choice of shape may vary greatly depending upon the desired results. Fluids 
exhibit differently in different shapes of microchannel. The flow rate and 
speed of the fluid in the microchannel are considered most meaningful for 
its suitability in any application. The chemistry of blood remained the same 
in all micro-channels [10]. 


Four micro-channels were subjected for the fluid optimization in the 
current study, which are briefly discussed below. 


1.1.1. Sinusoidal Microchannel. A sinusoidal shape microchannel 
comprised of periodic repetitions of oscillations [11]. Sinusoidal 
microchannels of defined wavelength and amplitude have been employed 
in previous research for heat dissipation and mass transformation [12]. 
Studies have showed that sinusoidal micro-channels exhibit excellent flow 
conditions for the blood-like fluid. Blood flow at curves tends to lose its 
streamlined flow; however, regained its normal tendon after the curves. 
Such micro-channels have been employed for biomedical implantations for 
varicose veins. Such micro-channels possessed excellent qualities for 
dissipation of heat because of the morphology [13]. Curvatures helped to 
resist and streamlining the flow by making it easier to homogenize in the 
fluids. 


Input 


Output 


Figure 1. Sinusoidal Microchannel 
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1.1.2. Spiral Microchannel. A coil-like microchannel called a spiral 
would have more complexity with increasing numbers of twists. Therefore, 
the current study’s application demanded minimum twist in channels. 
Continuous centrifugal force made the blood lose its streamline flow [14]. 
Morphology of retinal veins can easily afford the implantation of spiral 
micro-channels. Blood rheology easily allowed the use of a spiral channel 
for implantation purposes [15]. Such micro-channels have been widely 
employed in mass transit systems [16]. 


Input 


Figure 2. Spiral Microchannel 


1.1.3. U-Shaped Microchannels. Researchers have been employing U- 
Shaped microchannels in fluid dynamics effectively [17]. A big advantage 
of a U-shape is its ability to adjust its length as per need [18]. The single 
curve in such a channel tends to offer minimum parametric changes as 
compared to other shapes [19]. However, Its biomedical application is well 
investigated in previously done research [20]. 


Input Output 


Figure 3. U-Shaped Microchannel 
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1.1.4. Curvilinear Microchannel. A microchannel comprising of 
curving boundaries or edges is called a curvilinear microchannel. Blood 
particles in its plasma offer the same behaviour to curvilinear morphology 
fabricated microchannel as that of natural response to veins [21]. More 
prone to extendable flexibility, its suitability for biomedical implants is 
favourable [22]. With the increased number of curves, the cell separation 
task can be obtained from such micro-channels [23]. 


Input 


Output 


Figure 4. Curvilinear Microchannel 
1.2. Microchannels in Humans 


Biological systems in human organs, namely, brain, liver, lungs, kidney, 
heart, and eyes contain many micro-channels majorly performing the criti- 
cal functions of mass and heat transfer [24]. Cells can also migrate via mi- 
cro-channels in a body structure [25]. It is believed that a comprehensive 
microchannel mass transfer mechanism work in human biological systems. 
Human tissues contain a complex network of micro-channels, which can be 
used to regulate routine tissue operations [26]. Primarily, blood composed 
of two major components; plasma and red blood cells. Plasma is a Newto- 
nian fluid and red blood cells are intransigent spherical particles [27]. Blood 
bank viscosity is based on the percentage of blood that is made up of red 
blood cells. Characteristics of flowing blood like blood pressure, velocity, 
flow rate, the flow pattern is affected by the structure of the microchannel. 
Naturally occurring micro-channels in human biological system vary in 
morphology, which depends upon the function they are performing [28]. 
However, most micro-channels of human biological systems are elastic and 
can assume different morphology with the varying position of the human 
body [29]. Similarly, the diameter of this microchannel also responded to 
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the varying volume and pressure of blood and adjusted accordingly to tol- 
erate the instant conditions. 


1.3. Retinal Vein Occlusion 


The lens of the human eye focuses on all coming rays to form an image 
on the inner side of the eye sphere to a place called the retina [30]. Nerve 
cells surrounding the retina convert the coming light into signals, which are 
transmitted to the brain via optical nerve for image recognition. Blockage 
of blood flow in the retinal vein because of a clot is called retinal vein oc- 
clusion (RVO). RVO is one of the most commonly occurring retinal vascu- 
lar disease with an estimated 16 million effete globally [31]. A person af- 
fected by RVO may suddenly lose partial or full vision permanently without 
pain. So far iRVO is incurable; therefore, damage once done is irreversible 
but ophthalmologists can treat patients to stop further deterioration. Indi- 
viduals suffering from diabetes and high blood pressure are more prone to 
RVO [32]. Other commonly found reasons for RVO include hypertension, 
high cholesterol, smoking, glaucoma, and excessive weight [33]. 


Furthermore, RVO is additionally categorized into two types; Central 
Retinal Vein Occlusion (CRVO) and Branch Retinal Vein Occlusion 
(BRVO). In CRVO the main retinal vein carrying blood is blocked/raptured, 
whereas BRVO affects a small branch of the main retinal vein. Both central 
and branch retinal veins have diameters in the range of um with the central 
vein having a larger diameter as compared to the other branches. In certain 
cases, blocked vessels caused the leakage of jelly-like fluid into the retina 
resulting in small spots. These spots hampered the transmission of optical 
signals, which can be further used for the formation of an image. In other 
cases, some new unusual blood vessels may be caused by the retina; a pro- 
cess called neovascularization. These newly developed vessels leak the 
blood or fluid giving rise to RVO [34]. 


Different imaging tests are used to see the footsteps of RVO in the hu- 
man eye, which include fluorescein Angiography, Optical Coherence To- 
mography, and Fundus Autofluorescence [35]. In retinal veins of human’s 
maximum blood velocity and average flow rate as 54 mm/s and 2.04 nano- 
liters/ sec, respectively. Retinal veins carry blood away from the retina are 
quite similar to the microchannels. 
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1.3.1. Current Treatments. Current ophthalmology techniques to treat 
RVO include anti-vascular endothelial growth factor drug injections, retina 
opening massage,  clot-busting medication, steroid treatment, 
Panretinal photocoagulation, focal laser treatment, and hyperbaric oxygen 
therapy [36]. These treatments deal with the issues arising from RVO to 
prevent further damage. 


1.3.2. Challenges. There is no treatment or medication, which can 
completely unblock retinal veins [35]. Until now, RVO is incurable; 
ophthalmologists only manage the damage and prevent further deterioration. 
Vision problems, such as blurring or complete loss of sight in one eye, may 
appear unexpectedly. Complications from RVO include ocular edoema and 
haemorrhage. These conditions can cause blindness if not treated at the right 
time. 


1.4. Biomedical Implants 


Biomedical implants referred to placing natural or fabricated devices or 
tissues inside or on the surface of the human body to accomplish specific 
tasks. There may be different reasons for biomedical implants, which may 
include a missing body organ for supplemental support and monitoring 
body functions [37]. Hence, the implanted devices or tissues may be re- 
moved when they are no longer needed for the process. Primarily, opera- 
tional surgery was required for the implantation. Most commonly im- 
planted devices included implantable cardioverter defibrillators (ICDs), ar- 
tificial hips, heart pacemakers, breast implants, spine screws, bone support- 
ing rods, artificial spinal discs, intrauterine devices, metal supports for trau- 
matic fracture repair, artificial knees, coronary stents, ear tubes, and artifi- 
cial eye lenses. 


The interdisciplinary nature of biomedical implantation covering vari- 
ous fields like physics, chemistry, electronics, bioengineering, biotechnol- 
ogy, biophysics, pharmaceutics, and others in conjunction with micro and 
nanotechnology made this a frequent choice to fix chronic syndromes by 
many researchers [38]. 


1.4.1. Bioengineered Microchannels. In bypass surgery and dialysis, 
bioengineered veins provided a definite advantage to replace faulty natural 
veins with fabricated veins. Micro-channels for implantation to treat vari- 
cose veins problem were reported in the research [39]. Different materials 
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can opt for microchannel fabrication depending upon their application. Ini- 
tially, metal and ceramic-based micro-channels got the attention of previous 
researchers [40]. Furthermore, Polymer-based micro-channels got the atten- 
tion of researchers because of their biomedical compatibility. Especially the 
polymers with non-toxic and elastic properties gained an obvious advantage 
over metal-based micro-channels [41]. Polyvinylchloride, polypropylene, 
polyurethane, polyamide, polyethylenterephthalate, polyetherimide, and 
other synthetic polymers are among tried options for the fabrication of bio- 
engineered veins. Polydimethylsiloxane (PDMS) with its non-toxic proper- 
ties has gained the attention of researchers in last decade. Therefore, the 
current study investigated fluidic simulation and optimization of 
microchannels for Retinal Vein Occlusion (RVO). 


Centrai - 

Retinal Weim 
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Figure 5. Proposed Microchannel Implantation to Treat RVO 
2. FUZZY TECHNIQUE 


Artificial intelligence software enabled the researcher to deal with com- 
plex problems uncomplicatedly. Contrary to binary logic fuzzy logic pro- 
vided a more accurate picture of the system by considering all possibilities. 
Instead of overlooking at the implications of minutia, the fuzzy technique 
comprehensively enveloped all possible happenings. Fuzzy logic was ben- 
eficial in managerial sciences, business research, economics, scientific re- 
search, and various other engineering and social sciences fields. 


In MEMS technology mainly micro-size components are dealt with, 
where simulations of possible happening allow to save time, money, and 
energy for the process. Researchers have applied fuzzy techniques in the 
field of microfluidics to manage the fabrication and operation of devices 
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[42]. The fuzzy simulation depicted fluid flow properties and related out- 
comes well in advance. Morphology and geometry of micro-channels and 
other research employed in microfluidic devices were being predicted by 
fuzzy and ANSYS simulation for perfect outcomes [20, 43—49]. 


2.1. Fuzzy Simulation of Microchannels for RVO 


In the current study, four types of micro-channels architecture sinusoi- 
dal, U-shaped, spiral, and curvilinear were subjected to fuzzy parametric 
estimations to analyze the suitability for biomedical implantation to treat 
RVO. This microchannel architecture was selected to see the morphology 
of retinal veins. Moreover, the material used for fabrication can augment 
the evaluation of perfect-fit choices. 


For this purpose, four Fuzzy Logic Control (FLC) models one for each 
microchannel was prepared. Membership functions for each model was de- 
fined and rules were prepared by using “if and then” statements. 


3. RESULTS AND DISCUSSION 


The outcomes are discussed below for all three sinusoidal, curvilinear, 
and U-shaped microchannel. 


3.1. Sinusoidal Microchannel 


For sinusoidal microchannel, Fuzzy Inference System (FIS) variables 
and their ranges are shown in Table 1. 


Table 1. FIS Variables and their Ranges for Sinusoidal Microchannel 


Output FIS Variables with 


Input FIS Variables and Their Ranges Their Ranges 


Reynold 


Pressure Diameter 
Number 


Flow Rate Velocity 


10-100 KPa 100-1000 um 100-200 1-3 nL/s 40-60 mm/s 


The output critical for a retinal vein is flow rate and velocity. The range 
of Reynolds numbers was taken to take all possibilities into account. 


School of Science $ UMT— 57 


Volume 7 Issue 3, 2023 me 


Fluidic Simulation and Optimization... 


N 
a 


8 b N 


Figure 6. 3D Representation of Variations of Velocity and Flow rate in 
Sinusoidal microchannel 


Flow rate showed a great variance in response to varying pressure. Only 
for pressure near 85 KPa flow rate remained in the acceptable value near 
2nL/s. Response of flow rate to variations in Reynold’s number and diam- 
eter of the microchannel was spread over a wide range, showing variation 
in inputs for each unit change. This variation in the flow rate was a response 
to minor changes in all three input parameters, which made this the least 
compatible for the intended purpose. Velocity variations follow the same 
pattern of changes, which were followed in response for the input changes 
in the microchannel simulation. 


3.2. Spiral Microchannel 


Fore spiral microchannel Fuzzy Inference System (FIS) variables and 
their ranges are shown in Table 2. 


Table 2. FIS Variables and their Ranges for Spiral Micro-channels 


Output FIS Variables with 
Their Ranges 
Pressure Diameter Friction Flow Rate Velocity 
20-100 KPa 100-800 um 0-1 N 1-3 nL/s 40-60 mm/s 


Input FIS Variables and Their Ranges 
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Figure 7. 3D Representation of Variations of Velocity and Flow rate in 
Spiral Microchannel 


For variations in pressure, the flow rate varies to 2-2.75 nL/S. Only for 
40-80 KPa pressure ranges, the flow rate remained in acceptable ranges near 
2.1 nL/s. whereas in larger diameters the flow rate tends to increase con- 
stantly. Moreover, flow rate and velocity showed the same response of var- 
iations in the diameter of the microchannel. However, velocity remained 
constant near 5 mm/s for the pressure range of 55-10 KPa. 


3.3. U-Shaped Microchannel 


Fore U-Shaped microchannel Fuzzy Inference System (FIS) variables 
and their ranges are shown in Table 3. 


Table 3. FIS Variables and their Ranges for U-Shaped Micro-Channel 


Input FIS Variables and Their Output FIS Variables with Their 


Ranges Ranges 
Pressure Diameter Velocity Flow Rate 
10-100 KPa 100-1000 um 1-3 nL/s 40-60 mm/s 
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For low diameters up to 600 um, the flow rate maintains a constant value 
of 2nL/s, however, an increased flow rate was noticed in the range of 600- 
950 um that goes up to 2.7 nL/s. Contrarily, all pressures flow rate remained 
constant near 2nL/s. Minor changes in pressure and diameter made velocity 
vary to a great extent. 


Figure 8. 3D Representation of Variations of Velocity and Flow rate in U- 
Shaped Microchannel 


3.4. Curvilinear Microchannel 


Fore U-Shaped microchannel Fuzzy Inference System (FIS) variables 
and their ranges are shown in Table 4. 


Table 4. FIS Variables and their Ranges 


Input FIS Variables and Their Ranges uput ET Y ariamas ang 


Their Ranges 
Pressure Diameter Friction Flow Rate Velocity 
20-100 KPa 100-800um  0-1IN 1-3 nL/s 40-60 mm/s 


Curvilinear micro-channels possess the quality of U-Shape and sinusoi- 
dal micro-channels. Moreover, its morphological appearance resembles the 
natural retinal veins, which gives it a clear advantage over the other types. 


From the graphical representation it is quite obvious that for a diameter 
in the range of 400-600um its flow rate and velocity exhibit the ideal opti- 
mization. The most suitable situation in the curvilinear microchannel is that 
in the pressure range 50-100 KPa flow rate remains constant near 2 nL/s. 
Moreover, the flow rate also maintains a constant value near 2nL/s when 
the diameter of the channel ranges between 400-100 um. Velocity for any 
certain diameter maintains a constant value, which is ideal for the proposed 
application. 
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Figure 9. 3D Representation of Variations of Velocity and Flow rate in 
Curvilinear Microchanne 


3.5. Comparative Discussion 


A comparison of most optimum values was figured out through fuzzy 
parametric estimations, which is given below in Table 5. Both the flow rate 
and the velocity were compared with the Mamdani values to find the per- 
centage errors. Fuzzy rule views provided us with the ease to see all possible 
parameters for the optimal outputs. 


Although the percentage error in all four cases was quite negligible, the 
actual blood velocity and average flow rate values 54 mm/s and 2.04 nano- 
liters/ sec, respectively fall close in all simulated microchannels. 


UMT— 61 


School of Science Rony 
Volume 7 Issue 3, 2023 = 


Fluidic Simulation and Optimization... 


Table 5. Flow Rate and Velocity Results 


Microchannel Quantit FUZZY Mamdani Joage 
Type 7 y Value Value Error 
: : Flow Rate 2.18 nL/s 2.15 nL/s 1.4% 
Sinusoidal - 
Velocity 50.7 mm/s 50.5 mm/s 0.4% 
; Flow Rate 2.19 nL/s 2.16 nL/s 1.4% 
Spiral ; 
Velocity 51.6mm/s 51.9mm/s 0.6% 
Flow Rate 2.0 nL/s 2.3 nL/s 1.5% 
U-Shape ; 
Velocity 50 mm/s 50.2 mm/s 0.4% 
Flow Rate 2.06 nL/s 2.04 nL/s 0.1% 
Curvilinear ; 
Velocity 51.8 mm/s 51.5 mm/s 0.6% 


3.6. Conclusion 


The current study aimed to provide a fluidic stimulation by using Fuzzy 
parametric estimations technique, which provided optimal behaviours for 
all four micro-channels. The output optimization of four simulated micro- 
channels revealed that the curvilinear micro-channel best suited the implan- 
tation to treat RVO. Blood flow rate and velocity are the least divergent in 
a curvilinear microchannel. Although the simulated values of curvilinear 
micro-channels are a little lower than the practical value, this small defi- 
ciency is negligible in the context of associated flow rate. For the three sim- 
ulations, the combination of two output parameters remained divergent 
from the actual case requirements. 
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ABSTRACT 


Type-II diabetes is the most common type of diabetes, which has affected 
more than 465 million people globally and has become the ninth leading 
cause of mortality. The current study aims to determine risk factors 
associated with type-II diabetes such as BMI, cholesterol levels, physical 
activity, and smoking along with co-morbidities associated with this 
disease. The data was collected from the province of Punjab, Pakistan. For 
this purpose, a survey was performed for investigating the prevalence of 
type II diabetes and associated risk factors combined with it. A sample size 
of 265 patients was collected and interviewed through a questionnaire who 
were observed having type-II diabetes. A questionnaire was designed to 
record these patients’ responses. The questionnaire contained sub- 
categorizes, such as participants characteristics, lifestyle, and comorbidities 
associated with the disease. The findings indicated that the mean BMI was 
25.62 in the studied cohort. However, 78% of the patients had a family 
history of diabetes. Moreover, 76% of the participants were reported to be 
non-smokers and 75.67% of the participants reported atypical features 
associated with this disease. It was also observed that 13 % of the patients 
with low BMI (less than 18.5) had gastrointestinal diseases and 14.72% of 
the patients with a BMI greater than 18.5 but less than 25 had hypertension, 
and 8.30 % of the patients with more than 30 BMI also had hypertension. 
High BMI was found to be a major risk factor associated with type-II 
diabetes in this study population. The patients observed with high BMI were 
also more prone towards comorbidities associated with this disease. 
However, data suggested that most of the patients had familial type-II 
diabetes. 
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1. INTRODUCTION 


Non-communicable diseases (NCDs) are non-infectious and non- 
transmissible among people. These are chronic diseases arising due to 
genetic, environmental, and physiological factors. Globally, in 2010 34.5 
million deaths were caused by these NCDs [1]. It has been found that more 
than 60% of the global deaths occurred due to NCD in Korea, Japan, and 
China by following the static lifestyle, increased fat and protein intake in 
the diet are the major cause of increased mortality rate [2]. The mortality 
rate in Korea due to NCDs increased from 39.4% in 1983-56.0% in 2011 
[3]. In 2009, heart diseases, cerebrovascular disease, and neoplasm were 
identified with the prevailing causes of death resulting in more than 50% 
of the total amount of death [4]. For minimizing NCDs there must be a 
comprehensive understanding of the complete disease spectrum starting 
from health promotion and risk factor promotion to the rehabilitation. 
Among all, risk factor monitoring has been proven to be an efficient strategy 
in reducing NCDs [5]. The major cause of NCD includes the consumption 
of tobacco and alcohol, unhealthy dietary intake, and a sedentary lifestyle. 
In China, Japan, and Korea the transition in the diet pattern is the major 
cause of the rapid growth of NCD. They transformed their diet from the 
traditional vegetable pattern to an unhealthy diet, such as fast food and oily 
food. The protection and treatment of risk factors is proven to be more 
efficient in reducing mortality due to NCDs. Several public health centers 
target risk factor for preventing NCDs. Among China, Japan, and Korea 
some valuable achievements in risk factors have been noticed in the past 
few years. They have improved their dietary lifestyle [6]. 


In the Western countries the major cause of cardiovascular disease is 
high BMI (body mass index) according to previous epidemiological 
research [7], while in Asian countries low BMI is a major cause of 
cardiovascular death [8]. Prior studies have also indicated that people can 
prevent themselves from diabetes type II by having a healthy lifestyle 
routine. Different lifestyle factors have contributed to type I diabetes [9]. 
The risk factors for diabetes include obesity, weight gain, and physical 
inactivity, which can independently act as a risk factor [10]. 


According to the International Diabetes Federation (IDF), it was 
estimated that approximately 463 million of the adults worldwide (20-79 
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years of age) were suffering from type-II diabetes mellitus and by 2045 this 
number is expected to rise by 700 million. Furthermore, it was identified 
that 79% of the diabetic adults were living in low income countries. 
Moreover, 1 in every 2 individuals are undiagnosed diabetes patients. 
Globally, diabetes mellitus has caused over 4.2 million deaths and 374 
million individuals have been found to have an increased risk of developing 
mellitus type-II diabetes. 


Pakistan falls under the top 10 countries for absolute increasing rate due 
to diabetes prevalence, which reached 17.1% in 2019. A comparison with 
previously reported figures was drawn in this study, which suggested a 
148% increased percentage of this disease. Additionally, over 19 million 
adults in Pakistan were estimated with having diabetes, while 8.5 million 
out of these were undiagnosed [11]. This study focused on the prevalence 
of diabetes with associated risk factors and some atypical features (the 
features which are unusual and not the part of diagnostic procedures), which 
may or may not be linked to diabetes mellitus in Punjab, Pakistan. 


2. METHODOLOGY 


The current study deployed a cross-sectional survey to investigate the 
prevalence of NCDs and associated risk factors combined with the lifestyle 
and dietary controls of diabetic patients. The survey was conducted in 
Punjab province of Pakistan having a population of 220,892,340 as of 2020 
[12]. The sample size was calculated using the formula as described earlier 
[13]. 


A questionnaire consisting of 27 questions was designed for the survey, 
which was divided into three parts, namely comprising sub-sections as 
participants characteristics, lifestyle-related questions, and questions- 
related to comorbidities. Participants’ characteristics consisted of several 
questions related to the general information about the participants, such as 
their gender, age, weight, height, and family history. Participant’s weight 
was converted into pounds similarly participant’s height was recorded in 
feet was converted into inches to calculate the BMI of every participant 
[14]. 


A population of 265 participants were chosen having all type-II 
diabetes. The participants were made aware regarding the questions 
included in the survey, with each question being verbally explained to them 
by the conducting team. Furthermore, a written consent form was first filled 
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out by all the participants. Most of the lifestyle-related questions had ‘yes’ 
or ‘no’ answers. The analysis was compiled on Microsoft Excel with the 
descriptive analysis function and co-relations were analysed by using the 
AND function in MS Excel. 


3. RESULTS 
3.1. Participants Characteristics 


Among the selected participants 39.8% of the participants were males. 
Most of the participants belonged to the age group 50-59(25%), while the 
age group 20-29 contained the least of the participants (1.14%). Table 1 
summarizes the characteristics of the participants. The mean height in 
inches of the participants was 62.37 + S.D (standard deviation) of 9.7 
inches. Similarly, the mean weight in pounds was 145.12 + S.D of 2.54 Ibs 
and mean BMI was found to be 25.62 + S.D of 8. It was also observed that 
78% of the participants had a family history of diabetes disease. 


Table 1. Participants’ Characteristics of the Study (n= 265) 


Indices Number Percentage 
Gender 
Males 105 39.8% 
Females 159 60.2% 
Age 
<=19 51 19.32% 
20-29 3 1.14% 
30-39 16 6.06% 
40-49 50 18.94% 
50-59 66 25.00% 
60-69 45 17.05% 
70=> 33 12.50% 
BMI 
<18.5 64 24.81% 
18.5-24.9 65 25.19% 
25-30 57 22.09% 
>30 12 27.91% 
Family History 
Yes 207 18% 
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Indices Number Percentage 
No 58 22% 
(Mean BMI = 25.62 + 8) 
(Mean Weight in Pounds = 145.12 + 2.54 = 2.54) 
(Mean height in inches = 62.37 + 9.7) 


3.2. Lifestyle Behaviour 


The questionnaire targeted the participants who spent their holidays, 
while evaluating their diet and physical activity, which was significantly 
related to their disease. The exercise was among the foremost asked 
question in which 43.3% of the participants answered positively that they 
were actively engaged in exercise, while 28.35% answered “a few days a 
week”, whereas the same percentage answered “No”. When asked to choose 
between the work type, sedentary, mobile, and moderately active”, 43.68% 
of the participants chose sedentary, while only 17.24% chose mobile. (Table 
2) 


Another question of self-evaluation was considered to ask the 
participants whether they have checked their daily sugar intake or not. The 
results indicated that 75.67% of the participants chose Yes, which was a 
positive response. Participants were also asked about their dinner and sleep 
timings to calculate the interval between these activities. It was indicated 
thatl8% of the participants went to sleep in 2 hours, while 42% of the 
participants went to sleep right after 2 hours of having dinner, whereas 32% 
the participants were those who took 2-5 hours of sleep. Additionally, there 
were 3% of the participants who took more than 5 hours of sleep. However, 
2% of the participants had insomnia. 


When asked about smoking habits, 14% of the participants were 
smokers, while 76% were non-smokers. Additionally, 2% of the 
participants claimed that they quit and 7% said that they smoked On and 
off, while 1% of the participants answered Sometimes. 


Table 2. Lifestyle Behaviour of the Type-II Diabetic Patients (n= 265) 


Indices Number Percentage 
Exercise 
Daily 113 43.30% 
A few Days in a week 74 28.35% 
No 74 28.35% 
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Indices Number Percentage 
Work Type 
Sedentary 114 43.68% 
Moderately Active 102 39.08% 
Mobile 45 17.24% 
Check on Daily sugar intake 
Yes 199 75.67% 
No 64 24.33% 
Time interval between dinner and sleep 
<2 hours 47 18% 
Insomnia 5 2% 
2 hours 111 42% 
2 to 5 hours 86 32% 
More than 5 Hours 7 3% 
Smoking 
Yes 37 14% 
No 202 16% 
Sometimes 2 1% 
Quit 5 2% 
On and off 18 1% 
Check on Daily Sugar Level 
Yes 150 57% 
No 113 43% 


3.3. Atypical Features 


Participants were asked about atypical symptoms, which they might 
have that may or may not be linked to their diabetic health profile. These 
are given in Table III. However, these patients were asked about dizziness 
or weakness, around 199 (75.67%) participants agreed that they had these 
symptoms. Whereas 143(54.37%) of the participants reported the feeling 
pain or a burning sensation in their body. Additionally, 189(71.86%) 
reported numb feet and 53(20.15%) participants reported foot sores, while 
body ulcers were also reported by 58 (22.05%) of the participants. 
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Table 3. Atypical Features Associated with Type-II Diabetic Patients 


(n=265) 
Indices n (Percentages) n (Percentages) 
Atypical Features Yes No 
dizziness/weakness 199(75.67%) 64(24.33%) 
Feel of pain or Burning 143(54.37%) 120(45.63%) 


Numb feet 


189(71.86%) 


74(28.14%) 


foot Sores 53(20.15%) 210(79.85%) 
Body ulcers 58 (22.05%) 205(77.95%) 
Amputation 29 (11.03%) 233(88.59%) 


3.4. Co-analysis 


Table 4 highlights the relationship between different categories of 
questions. A percentage of 6.79% of the participants were mentioned having 
all the symptoms as given in Table 3, while 80% of the participants had at 
least one of these symptoms (Table 3). When a link between lifestyle and 
work type of the participants was made, it was found that 16.98% of the 
participants had a sedentary lifestyle and did not exercise at all. The 
participants, which had a mobile lifestyle and exercised daily were observed 
to be 18(6.79%), while 19(7.17%) participants had a moderately active 
lifestyle and exercised a few days in a week. It was observed that patients 
with low BMI (less than 18.5) were indulged in stomach, liver, and intestine 
diseases with a percentage of 13%. Participants with BMI greater than 18.5 
but less than 25 had Hypertension (14.72%) and patients with more than 30 
BMI had hypertension as well (8.30%). 


Table 4. Co-analysis of Different Risk Factors in Type-II Diabetic Patients 
(n=265) 


Serial dices Participants 
No. (Percentage) 
1. Participants that have at least one of the 212(80.00%) 
mentioned symptoms 
2 Participants that have a sedentary lifestyle and do 45(16.98%) 
not exercise 
3. Participants that have a moderately active 19(7.17%) 


lifestyle and exercise a few days in a week 
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Serial indices Participants 
No. (Percentage) 
4. Participants that have a mobile lifestyle and 18(6.79%) 
exercise daily 
5. Participants that have all the mentioned 18(6.79%) 
symptoms 
6. Participants with <18.5 BMI and have stomach, 35(13%) 
liver or intestine diseases 
ds People with <25 BMI and have Hypertension 39(14.72%) 
8. People with >30 BMI and have hypertension 22(8.30%) 


Table 5 gives the prevalence of co-morbidities with diabetes. 
Hypertension was the most common disease with 81.89% of the participants 
suffering from it. Blood cholesterol was the second highest with 81.51% of 
the participants. Additionally, 58.49% of the participants reported heart 
disease. However, other (40.75%) reported diseases related to Stomach, 
intestine, and liver, while in a separate category of GI diseases, 40.75% of 
the participants were identified with these diseases. Among all only 2.26% 
of the participants reported having paralysis issue. 


Table 5. Prevalence of Co-morbidities with Type —II Diabetes 


HT 217 (81.89%) 
Cholesterol 216(81.51%) 
Heart Diseases 155(58.49%) 
Stomach Liver Intestine 80(30.19%) 
GI 108(40.75%) 
Paralysis 6(2.26%) 


4. DISCUSSION 


Non communicable diseases (NCDs) have been linked with effective 
lifestyle behaviours. It has been reported that a successful change in lifestyle 
caused a 58% reduction in the diabetic incidence rate [15]. However, in 
developing countries, more than 80% of the deaths have been reported as a 
result of cardiovascular diseases and diabetes [16]. Behavioural risk factors 
and physiological risk factors have been reported to be the primary causes 
for NCDs [17]. Physical activity has been inversely related to different 
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diseases, such as diabetes [18]. Nevertheless, even today, there are quite 
less solutions that deal with physical inactivity. Reportedly, data collected 
from Brazil reported that 61% of the adults did not reach satisfactory levels 
of physical activity [19]. 


Hayes, L., et al studied the prevalence of diabetes in Pakistan before and 
after the physical activity of Pakistani population and its risk factors, which 
are associated with diabetes. He indicated that Pakistanis are less physically 
active than European people [20]. This corroborates with the current study 
in which 7.17% participants were found to have a moderately active 
lifestyle and people who were indulged in exercising were more active. A 
good percentage of participants approximately 43.68% had a sedentary 
work type and 16.98% out of these did not considered physical exercising 
as a part of active lifestyle behaviour. 


Diabetes is also linked with sleep disturbances [21], which was actively 
found in this study, 2% of the participants reported insomnia, while 3% of 
the participants took more than 5 hours to sleep after taking dinner. It is 
worth mentioning that a considerable amount, approximately 32% took 2-5 
hours of sleep, which could either be a lifestyle choice or can also be the 
cause of sleep disturbances [22]. 5% of the participants had sleep disorders, 
which can also be linked to diabetes. 


Reportedly, this study identified that patients having diabetes were 
diagnosed with body ulcers and foot sores, which were 22.05% and 20.15%, 
respectively. According to some previous studies, there are serious 
complications among patients suffering from diabetes contributing to 
morbidity and mortality [23]. Noticeably, patients who were often 
prescribed inappropriate antibiotics worsened the situation leading to 
adverse effects, unlikely to heal. Diabetics with active ulcers have healing 
rates up to 75%. The amputation rate in this study’s findings is 11.03%, 
while 15-20% of the amputations are of prior studies [24]. 


An extensive number of the participants, approximately 71.86% 
reported numbness of feet, which can be an indication of diabetic 
neuropathy [25]. Furthermore, it was observed in this study that diabetic 
patients who are overweight (25-30 BMI) and obese (>30 BMI) suffered 
from hypertension. This association corresponds with the results of previous 
findings. Moreover, higher levels of BMI was linked to the increasing blood 
pressure of patients and the gradient observed strongly suggested that there 
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could be a direct effect of BMI on blood pressure disregarding other clinical 
factors [26]. The results of these findings indicated a mean BMI of 25.62 + 
8. with males having a mean BMI of 25.69 and females having a BMI of 
29.58 as compared to the body mass index 24.1 and 24.7, respectively of 
diabetic patients in Japan [27]. 


Several studies conducted in Asia identified Stroke and coronary heart 
diseases are also associated with comorbid conditions of diabetes. 
Moreover, these studies indicated diabetes as a leading cause of death in 
patients and with reported stroke rate of 42%, respectively [28]. However, 
no valuable association was found between BMI and stroke in the 
epidemiologic report of Japan [29]. Smoking has been reported as a 
significant associated risk factor for developing diabetes with 44% of the 
increased risk rate [30]. Cigarette smoking induces insulin resistance and 
inadequate compensatory insulin secretion responses. 


4.1. Conclusion 


The study was conducted on a population sample of 265 participants, 
which were suffering from the diabetes. Noticeably, high BMI was found 
to be a major risk factor associated with type-II diabetes in this study 
population. Moreover, patients having a high BMI were also more prone 
towards comorbidities associated with this disease. Markedly, proper diet 
and exercise can reduce the risk of type-II diabetes among patients. 
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ABSTRACT 


Actinomycetes abundantly present in soil are renowned for various 
secondary metabolites production with a wide range of use in diverse fields 
of life. Many reports are available on the use of actinomycetes for the 
regulation of pathogenic microbes in plants. For this purpose, this study 
used thirteen Actinomycetes strains, which were isolated from the fields of 
Piplan (district Mianwali) Punjab, Pakistan. However, out of thirteen, only 
seven were found active for the production of both primary and secondary 
screening. The active seven actinomycetes (AWI, AW2, AW3, AW4, 
AM1, AC1, and AC2) were separated for their microbial potential against 
pathogenic strains (Escherichia coli and Staphylococcus aureus) isolated 
from the same agriculture soil. About 42.85% of isolated actinomycetes 
(three out of seven) showed antagonistic properties against S. aureu and E. 
coli in primary screening. Thus, it was noticed that the three strains (AW3, 
AM1, and AC1) had great antimicrobial potential, which showed the most 
promising results. Whereas in Secondary screening only ACI showed the 
best result with ethyl acetate extract as compared to ethanol, methanol, and 
chloroform. Hence, It was concluded that actinomycetes isolates had great 
potential for antibacterial activity and they can be used in agriculture 
industry for the regulation of plant pathogens. 


Keywords: Actinomycetes, antimicrobial, metabolites, piplan, primary 
screening, secondary screening 


1. INTRODUCTION 


Actinomycetes belong to a group of bacteria placed in order in 
Actinomycetales, having a distinguished taxonomic group in the bacterial 
domain; however, with fungi, they share certain morphological features [1]. 
They are aerobic, spore formers, and Gram-positive. In soil, they create 
thread-like filaments [2]. They display filamentous growth, producing 
substrate or aerial mycelium, occasionally coccoid to rod in shape [3]. 
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Actinomycetes DNA has high G+C content in the range of 57-75% [4]. Out 
of eighteen main lineages of the bacterial domain, Actinomycetes are one 
of the prime taxonomic units [5]. They cultivate as hyphae and are 
characteristically capable of the production of soil earthy smell [4, 6]. 
According to [7] around 100 actinomycetes genera are found in soil and 
several species, which are primarily found as soil inhabitants. 


However, reportedly in the previous literature, Actinomycetes existed 
as an appearing group of microorganisms, which were extensively dispersed 
in natural ecosystems around the world [8, 9]. Furthermore, the species had 
been found in a divergent aquatic environment containing deep ocean 
sediments [10—12]. 


It is also known that the Actinomycetes are producers of 
many secondary metabolites [13]. In actinomycetes, about 70-80% of 
secondary metabolites are produced by Streptomyces genera (almost 500 
species); while other genera have less contribution, for example, 
Micromonospora, Saccharopolyspora, Actinoplanes, and Amycolatopsis 
[14]. They have a broad range of utilization within the agricultural, medical, 
and pharmaceutical industries. In medicine, secondary products are used 
as antibiotics, antitumors, antiviral, and anti-infection agents [15-17]. 
Furthermore, actinomycetes would manufacture plant growth promoters to 
boost plant growth and facilitate the institution of plants under stressed 
conditions [18, 19]. Secondary metabolites produced by actinomycetes 
facilitate other microbes by forming mutualistic and symbiotic interactions 
with them to protect their host from pathogenic microbes [20]. There is 
scarce data from fields of Piplan (district Mianwali) about the isolation of 
ascomycetes and their antimicrobial activity from the rhizospheric region 
of the crops. Therefore, the current study aimed to investigate the isolation 
of actinomycetes from (fields of piplan) and the use of these strains against 
pathogenic strains, which were isolated from crops such as maize (Zea 
mays), cotton (Gossypium), and wheat (Triticum aestivum). The current 
study would significantly highlight the applications of isolated 
actinomycetes for the control of plant pathogens strains. 
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2. MATERIALS AND METHODS 
2.1. Collection of Soil Samples 


Soil samples were collected from the rhizospheric region of the crops 
such as wheat (Triticumaestivum), cotton (Gossypium), and maize (Zea 
mays) for the isolation of actinobacteria. 


The soil samples were sited in polythene bags (sterile) during the month 
of November 2017. 


The samples were desiccated at ambient temperature for continuous 2 
days and sieved. The sieved soils were then used for the actinomycetes 
isolation. 


2.2. Isolation and Characterization of Actinobacteria 


For isolation purpose, enrichment of Actinomycetes in the samples was 
carried out by physical [21] and chemical treatment [22]. Following 
physical and chemical treatment 1 gm of soil sample was serially diluted 
and spread on ISP media supplemented with 50 g/ml antifungal agent 
nystatin [23]. The isolated actinobacteria (13 strains) were sub-cultured and 
incubated for 7—12 days at 28°C. However, only seven selected colonies of 
actinobacteria were noticeably characterized both morphologically and 
biochemically [24]. 


2.3. Test Organisms 


Pathogenic soil bacterial strains, such as S. aureus and E. coli were 
isolated from the same soil samples, which were used to assess the 
antagonistic activity of actinomycetes. Glycerol stock was prepared to store 
bacterial strains. It can be used by subculturing when required. 


2.4. Primary Screening for the Antibacterial Activity of Actinomycetes 


In primary screening, a modified cross-streak process was used to assess 
antimicrobial activity [25]. From each plate, an isolated colony of 
actinomycetes was picked and streaked on the modified nutrient agar 
surface. 


A single streak of actinomycetes isolates was created for each isolate, 
which passed down the middle of the plate and incubated at 28°C for 5-7 
days [26]. An upright ribbon growth of the actinomycetes was observed on 
the plates, followed by the inoculation within 24 hours old culture 
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pathogenic soil bacterial strains i.e. E. coli and S. aureus. Then, the plates 
were further incubated for 24 h at 37 °C [27]. The inhibition zone was 
measured and recorded after the incubation. 


2.5. Secondary Screening for the Antibacterial activity of 
Actinomycetes 


The Agar well diffusion process along with crude extracts of ethanol, 
ethyl acetate, chloroform, and methanol was used in the secondary 
screening of the isolate. Submerged fermentation bioactive compound 
production was done by Chaudhary et al. [2]. A conical 250 ml flask was 
taken with 50 ml of ISP broth, which was inoculated with Actinomycetes 
isolates. After the incubation was done at 30°C in a shaking incubator at 
150-rpm rotation for continuous 7 days. When fermentation was done and 
centrifuged at medium twice for 10 min at 10,000 rpm to take out cells and 
debris and harvested them for fermented broth. Collected crude extracts 
were tested against test microorganisms (S. aureus and E.coli) via a method 
of agar well diffusion [28, 29]. Both test microorganism and cell 
concentrations were adjusted at 0.5 McFarland turbidity standards. A 
sterilized cotton swabs were used to inoculate the strains on nutrient agar 
plates. A sterilized micro tip (1000 ul) was used to bore wells in the plates. 
Then, the wells were poured with 20, 50, and 100 ul of each crude extract. 
The plates were incubated at 37°C for 24 h [2]. After incubation, the zone 
of inhibition was measured and noted for the record. 


3. RESULTS AND DISCUSSION 


In the last few years, actinomycetes seemed extensively studied in 
various undiscovered locations in distinct parts of the globe (including 
Pakistan). However, there are no reports of isolation of actinomycetes from, 
District: Mianwali Punjab (Pakistan). Therefore, attempts have been made 
to isolate actinomycetes from this undiscovered area to find new species. 
From Soil samples, 7 actinomycetes strains (out of 13) were isolated and 
labeled as Actinomycete wheat 1 (AW1), Actinomycete wheat 2 (AW2), 
Actinomycete wheat 3 (AW3), Actinomycete wheat 4 (AW4), 
Actinomycete maize 1 (AM1), Actinomycete cotton 1(AC1), and 
Actinomycete cotton 2 (AC2). In the current study overall, 53.8 % isolates 
of actinomycetes were found to be gram-positive organisms and the results 
were similar with Subhan et al. [30] who identified 55 % isolates of 
actinomycetes as gram-positive. Furthermore, two pathogenic bacterial 
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strains (E.coli and S. aureus) were also isolated, which were characterized 
both morphologically and biochemically (see Table 1). 


Table 1. Morphological and Biochemical Characteristics of Bacterial 
Strains Isolated from Agricultural Soil 


Sr.No. GS S M C Sp F C Orn MR VP Cit Nt U HS 
S. aureus + Cocci - - -_  - +- - + + + + + - 
E. Coli - Rods + - - + + - + + - = go os - 


Note. GS: Gram Staining; S: Shape; M: Motility; C: Capsule; Sp: Spore; F: 
Flagella; C: Catalase; O: Oxidase; In: Indole Production; MR: Methyl Red; 
VP: Voges-Proskauer; C: Citrate Utilization; Nit: Nitrate Reduction; U: 
Urease; H2S: H2S Production 


The morphological characterization of all the seven isolates showed dry 
and rough to the powdery and smooth texture of isolates. Moreover, the 
colonies were chalky white, dry, nodular, and sticky on isolation agar of 
actinomycetes. The morphological and biochemical characteristics 
suggested that the isolates were representative of genus Actinomycetes (see 
Table 2). The characteristics of the genus were also in line like the previous 
findings of Gurung et al. [31] and Anwar et al. [32]. 


Table 2. Biochemical Characteristics of Actinomycetes Bacterial Strains 


Isolated strains of Actinomycetes 
AWI1 AW2 AW3 AW4 ACI AC2 AMI 


Biochemical Test 


Catalase + - + + + + + 
Citrate Utilization + - + + + + 2 
Urease + + - 3 + z + 
Methyl Red + + + + + + + 
Voges-Proskauer - + + + = Š E 
Melanin Production + + + - + - + 
H2S Production - - + : 2 2 z 
Nitrate Reduction + + + + + - + 
Motility - - : = = = : 
Starch Hydrolysis + + + + + + + 


Note. AW1: Actinomycete wheat 1; AW2: Actinomycete wheat 2; AW3: 
Actinomycete wheat 3; AW4: Actinomycete wheat 4; AM1: Actinomycete 
maize 1; AC1: Actinomycete cotton 1 and AC2: Actinomycete cotton 2 


In the current study, a soil sample was collected from the surface and 
rhizospheric region of crops at the depth and it has been found that as a 
result of appropriate pH and water condition in depth, actinomycetes 
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quantity at depth is quite higher than the surface of the soil and similar 
findings were also reported Basavaraj et al. [33]. 


For antimicrobial potential, strains of actinomycetes were assessed. By 
perpendicular streak method, primary screenings were performed [29]. 
Firstly, the perpendicular streak method (primary screening) was used to 
test actinomycetes strains and it was concluded that three (AW3, AM1, and 
AC1) out of seven strains (AW1, AW2, AW3, AM1, AC1, AW4, and AC2) 
(see Table3) showed activity against test S. aureus and E. coli (gram- 
positive and gram-negative bacteria). This broad-spectrum activity of 
actinomycetes may be attributable to the production of more than one 
antimicrobial compound, which compels the strains more effectively 
against the both tested bacteria [31]. An overall 42.5% activity of 
actinomycetes (three out of seven strains) was recorded and founding was 
in close agreement with Seipke et al. [34] who stated 45 % of isolates; 
whereas Remya and Kumar [35] reported that 47 % of isolates showed 
broad-spectrum activity in the primary screening. However, a lower 
antimicrobial activity of actinomycetes of 38 % (51 out of 134 
actinomycetes) were compared to this study, which was reported by 
Sharma et al. [36]. Similarly, findings of low antimicrobial activity of 
actinomycetes were reported by [37, 38]. Moreover, Belyagoubi et al. [39] 
reported a very high level of 72.86 % of antimicrobial activity of 
actinomycetes as compared to this study’s findings. 


Table 3. Primary Screening against Bacterial Isolates 


Actinomycetes and zone of inhibition (mm) 
AW1 AW2 AW3 AW4 AMI ACI AC2 
Staphylococcus aureus 10.5 0 10.5 0 20.5 20.6 0 
Escherichia coli 0 10.5 20.5 0 10.5 20.8 0 


Bacterial Isolate 


Table 4. Secondary Screening Results 


: Actinomycetes 
Bacterial Isolate AW} AMI ACI 
Staphylococcus aureus z = + 
Escherichia coli - = $ 


For secondary screening, three actinomycetes (AW3, AM1, and AC1) 
isolates were selected based on the primary screening (see Table 4). The 
Agar well method was used to perform secondary screening of active 
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isolates [28, 29]. A finding of secondary screening was conclusive. In all 
three actinomycete isolates only Actinomycetes ACI showed great 
potential for antimicrobial activity as compared to those strains, which 
exhibit activity in primary screening but did not give antimicrobial activity 
in secondary screening. The reason behind this is maybe when growing on 
a solid medium the modification in the morphology of actinomycetes 
(filamentous mycelia) while in liquid broth (fragmenting mycelia) [28] and 
mostly actinomycetes are poor fermenters. Actinomycetes release active 
compounds, which became inactive or stick to the liquid medium 
component or in broth they modified chemically is also the possibility [31, 
40]. 


Table 5. Zone of Inhibition in mm with Ethyl Acetate Extract 


Zone of inhibition (mm) 


Bacterial Isolate 


20uL 50 uL 100 uL 
Staphylococcus aureus 3.1 5.5 12 
Escherichia coli 2.9 5l 10 


Different solvents were used for the isolation of antimicrobial 
metabolites including ethyl acetate, ethanol, chloroform, and methanol from 
the broth. This method is called a solvent extraction method. Only ethyl 
acetate extract of actinomycetes AC1 exhibited antimicrobial potential 
against Staphylococcus aureus and Escherichia coli. The recorded activity 
of AC1 against Staphylococcus aureus was measured by inhibition zone as 
3.1, 5.5 and 12 mm with 20, 50, and 100 mg/L concentrated extracts of ethyl 
acetate respectively. However, the recorded activity of AC1 against E.coli 
was measured by inhibition zone as 2.9, 5.1, and 10 mm with 20, 50, and 
100 mg/L concentrated extracts of ethyl acetate respectively (Table:5). It 
was noticed in the present study, that fermented broth with other solvents 
did not extract any antimicrobial metabolite as no inhibition zone was 
observed with any concentration. Gurung et al. [31] findings with other 
solvents suggest the failure of metabolite extraction may be due to the 
existence of polar functional group in a secondary metabolite that form 
metabolite soluble in water and insoluble in solvent, the use of inappropriate 
solvents and inadequate shaking of the mixture. 


3.1. Conclusion 


The current study showed that isolates of Actinomycetes, found in 
Piplan soil of Mianwali, had the potential to use it as sources for the novel 
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antibacterial compounds against pathogenic microorganisms of soil. Isolate 
ACI exhibited the highest activity in secondary screening against both 
gram-positive and gram-negative test bacteria. However, further work in 
this research is required about AC1 for its practical application in this field. 


REFERENCES 


1. 


Hirsch CF, Christensen DL. Novel method for selective isolation of 
actinomycetes. Appl Environ Microbiol. 1983;46(4):925—929. 
https://doi.org/10.1128/aem.46.4.925-929.1983. 


Chaudhary HS, Yadav J, Shrivastava AR, Singh S, Singh AK, Gopalan 
N. Antibacterial activity of actinomycetes isolated from different soil 
samples of Sheopur (A city of central India). J Adv Pharma Technol 
Res. 2013;4(2):118—123. _ https://doi.org/10.4103%2F2231-4040. 
111528 


. Adegboye MF, Babalola OO. Taxonomy and ecology of antibiotic 


producing actinomycetes. Afr J Agric Res. 2012;7(15):2255—22561. 
https://doi.org/10.5897/AJARX11.071 


Boroujeni ME, Arijit D, Prashanthi K, Sandeep S, Sourav B. Enzymatic 
screening and random amplified polymorphic DNA fingerprinting of 
soil streptomycetes isolated from Wayanad district in Kerala, India. J 
Biol Sci. 2012;12(1):43—50. https://doi.org/10.3923/jbs.2012.43.50 


. Jeffrey LS. Isolation, characterization and identification of 


actinomycetes from agriculture soils at Semongok, Sarawak. Afr J 
Biotechnol. 2008;7(20): 3697—3702 


Sprusansky O, Stirrett K, Skinner D, Denoya C, Westpheling J. The 
bkdR gene of Streptomyces coelicolor is required for morphogenesis 
and antibiotic production and encodes a transcriptional regulator of a 
branched-chain amino acid dehydrogenase complex. J Bacteriol. 
2005; 187(2):664—67 1. https://doi.org/10.1128/jb.187.2.664-67 1.2005 


Yamamoto Y, Kouchiwa T, Hodoki Y, Hotta K, Uchida H, Harada KI. 
Distribution and identification of actinomycetes lysing cyanobacteria in 
a eutrophic lake. J Appl Phycol. 1998;10:391—397. https://doi. 
org/10.1023/A: 1008077414808 


Scientific Inquiry and Review 


»— SIR 


Volume 7 Issue 3, 2023 


Zamir et al. 


10. 


11. 


12. 


13. 


14. 


15. 


16. 


Singh R, Dubey AK. Diversity and applications of endophytic 
actinobacteria of plants in special and other ecological niches. Front 
Microbiol. 2018;9:e1767. https://doi.org/10.3389/fmicb.2018.01767 


Tandale A, Khandagale M, Palaskar R, Kulkarni S. Isolation of pigment 
producing actinomycetes from soil and screening their antibacterial 
activities against different microbial isolates. Int J Curr Res Life Sci. 
2018;7(6):2397—2402. 


Arifuzzaman M, Khatun MR, Rahman H. Isolation and screening of 
actinomycetes from Sundarbans soil for antibacterial activity. Afr J 
Biotechnol. 2010;9(29):4615—4619. 


Ayoubi H, Mouslim A, Moujabbir S, et al. Isolation and phenotypic 
characterization of actinomycetes from Rabat neighborhood soil and 
their potential to produce bioactive compounds. Afr J Microbiol Res. 
2018;12(8):186—191. https://doi.org/10.5897/AJMR2017.8761 


Kirby BM, Roes-Hill ML, Cary SC, Burton SG, Tuffin IM, Cowan DA. 
Actinobacterial Diversity Associated with Antarctic Dry Valley Mineral 
Soils. In: de Bruijn FJ, ed. Handbook of Molecular Microbial Ecology 
IT; Metagenomics in Different Habitats. Wiley-Blackwell; 2011:125— 
33. https://doi.org/10.1002/9781118010549.ch13 

George M, Anjumol A, George G, Hatha AM. Distribution and 
bioactive potential of soil actinomycetes from different ecological 
habitats. Afr J Microbiol Res. 2012;6(10):2265—2271. https://doi.org/ 
10.5897/AJMR11.856 


Berdy J. Bioactive microbial metabolites. J Antibiot. 2005;58(1):1—26. 
https://doi.org/10.1038/ja.2005.1 

Chen Y, Zhou D, Qi D, Gao Z, Xie J, Luo Y. Growth promotion and 
disease suppression ability of a Streptomyces sp. CB-75 from banana 


rhizosphere soil. Front Microbiol. 2018;8:e2704. https://doi.org 
/10.3389/fmicb.2017.02704 


Gallagher KA, Fenical W, Jensen PR. Hybrid isoprenoid secondary 
metabolite production in terrestrial and marine actinomycetes. Curr 
Opinion Biotechnol. 2010;21(6):794—800. https://doi.org/10.1016/). 


copbio.2010.09.010 


UMT—- 91 


School of Science D 
Volume 7 Issue 3, 2023 wwe 


Screening of Actinomycetes Isolated... 


17. 


18. 


19. 


20. 


21. 


22. 


23; 


24. 


25. 


26. 


27. 


Igarashi Y. Screening of novel bioactive compounds from plant- 
associated actinomycetes. Actinomycetologica. 2004;18(2):63-66. 
https://doi.org/10.3209/saj.18_63 


Bonjar GS. Broadspectrim, a Novel Antibacterial from Streptomyces 
sp. Biotechnology. 2004;3:126-130. 


Srivastava S, Patel JS, Singh HB, Sinha A, Sarma BK. Streptomyces 
rochei SM 3 induces stress tolerance in chickpea against Sclerotinia 
sclerotiorum and NaCl. J Phytopathol. 2015;163(7-8):583-592. 
https://doi.org/10.1111/jph.12358 


Seipke RF, Kaltenpoth M, Hutchings MI. Streptomyces as symbionts: 
An emerging and widespread theme? FEMS Microbiol Rev. 
2012;36(4):862-876. https://doi.org/10.1111/).1574-6976.2011. 
00313.x 


Seong CN, Choi JH, Baik KS. An improved selective isolation of rare 
actinomycetes from forest soil. J Microbiol. 2001;39(1):17-23. 


Oskay M. Comparison of Streptomyces diversity between agricultural 
and non-agricultural soils by using various culture media. Sci Res 
Essays. 2009;4(10):997—1005. 


Taechowisan T, Peberdy JF, Lumyong S. Chitinase production by 
endophytic Streptomyces aureofaciens CMUAc130 and its antagonism 
against phytopathogenic fungi. Ann Microbiol. 2003;53(4):447-462. 


Ali B, Sabri AN, Ljung K, Hasnain S. Quantification of indole-3-acetic 
acid from plant associated Bacillus spp. and their phytostimulatory 
effect on Vigna radiata (L.) World J Microbiol Biotechnol. 


2009;25:519-526. https://doi.org/10.1007/s11274-008-99 18-9 


Holt JG. Bergey’s manual of systemic bacteriology. UK: Cambridge 
University Press; 1989. 


Hemashenpagam N. Purification of secondary metabolites from soil 
actinomycetes. Int J Microbiol Res. 2011;3(3):148-156. 


Bizuye A, Moges F, Andualem B. Isolation and screening of antibiotic 
producing actinomycetes from soils in Gondar town, North West 
Ethiopia. Asian Pacific J Trop Dis. 2013;3(5):375—381. https://doi.org 
/10.1016/S2222-1808(13)60087-0 


Scientific Inquiry and Review 


2—SIR 


Volume 7 Issue 3, 2023 


Zamir et al. 


28. 


29. 


30. 


31. 


32. 


33. 


34. 


35. 


36. 


Budhathoki S, Shrestha A. Screening of actinomycetes from soil for 
antibacterial activity. Nepal J Biotechnol. 2020;8(3):102—110. 
https://doi.org/10.3126/njb.v813.33664 


Wahab A, Shumaila S, Subhan SA, Ali ST, Mujahid TY. Isolation and 
identification of actinomycetes isolated from Karachi soil and screening 
of antimicrobial compounds. Int J Curr Res. 2015;7(2):12760-12765. 


Subhan SA, Wahab A, Mujahid TY, Abbas T, Khan I, Idrees S. 
Screening of actinomycetes from indigenous soil for production of 
extracellular metabolites. Int J Curr Res. 2015;7:12078—12083. 
https://doi.org/10.1016/S2222-1808(13)60087-0 


Gurung TD, Sherpa C, Agrawal VP, Lekhak B. Isolation and 
characterization of antibacterial actinomycetes from soil samples of 
Kalapatthar, Mount Everest Region. Nepal J Sci Technol. 2009; 10:173- 
182. https://doi.org/10.3126/njst.v1010.2957 


Anwar S, Ali B, Sajid I. Screening of rhizospheric actinomycetes for 
various in-vitro and in-vivo plant growth promoting (PGP) traits and for 
agroactive compounds. Front Microbiol. 2016;7:e1334. 
https://doi.org/10.3389/fmicb.2016.01334 


Basavaraj KN, Chandrashekhara S, Shamarez AM, Goudanavar PS, 
Manvi FV. Isolation and morphological characterization of antibiotic 
producing actinomycetes. Trop J Pharma Res. 2010;9(3):231—236. 
https://doi.org/10.4314/tjpr.v913.56282 


Seipke RF, Kaltenpoth M, Hutchings MI. Streptomyces as symbionts: 
an emerging and widespread theme? FEMS Microbiol Rev. 


2012;36(4):862-876. https://doi.org/10.1111/).1574-6976.2011. 
00313.x 
Remya M, Vijayakumar R. Isolation and characterization of marine 


antagonistic actinomycetes from west coast of India. Med Biol. 
2008;15(1):13-19. 

Sharma D, Kaur T, Chadha BS, Manhas RK. Antimicrobial activity of 
actinomycetes against multidrug resistant Staphylococcus aureus, E. 
coli and various other pathogens. Trop J Pharma Res. 2011;10(6):801— 
808. https://doi.org/10.4314/tjpr.v 1016.14 


UMT— 93 


School of Science Rost 
Volume 7 Issue 3, 2023 wwe 


Screening of Actinomycetes Isolated... 


37. 


38. 


39. 


40. 


Pandey B, Ghimire P, Agrawal VP. Studies on the antibacterial activity 
of the Actinomycetes isolated from the Khumbu region of Nepal. J Biol 
Sci. 2004;23:44-53. 


Rana S, Salam MD. Antimicrobial potential of actinomycetes isolated 
from soil samples of Punjab. India. J Microbiol Exp. 2014;1(2): 63—68. 
https://doi.org/10.15406/jmen.2014.01.00010 


Belyagoubi L, Belyagoubi-Benhammou N, Jurado V, et al. 
Antimicrobial activities of culturable microorganisms (actinomycetes 
and fungi) isolated from Chaabe Cave, Algeria. Int J Speleol. 
2018;47(2):189-199. https://doi.org/10.5038/1827-806X.47.2.2148 


Aedms H, Sheir DH, Eldewany AI. Production of antimicrobial agent 
from marine bacteria Isolated from Mediterranean. Aust J Basic Appl 
Sci. 2011;5(5):121-128. 


Scientific Inquiry and Review 


4— SIR 


Volume 7 Issue 3, 2023 


Scientific Inquiry and Review (SIR) 
Volume 7 Issue3, 2023 
ISSN p: 2521-2427, ISSN œ: 2521-2435 


Homepage: https://journals.umt.edu.pk/index.php/SIR 


Title: 


Author (s): 


Affiliation (s): 


DOI: 
History: 


Citation: 


Copyright: 


Licensing: 


Conflict of 
Interest: 


Synthesis, Properties, and Applications of Carbon Nanotubes: An 
Overview 


Shabbir Hussain!, Uzma Akbar?, Muhammad Ahmad?, Muhammad Ibrar*, Zulfiqar 
Ali, Muhammad Waqas‘, Sheikh Asrar Ahmad®, Syed Mustansar Abbas’, Habib 
Ullah*, Sehrish Anwar* 


'Khwaja Fareed University of Engineering and Information Technology Rahim Yar Khan, 
Pakistan 

*Minhaj University Lahore, Pakistan 

3University of Education, Lahore, Pakistan 

‘Lahore Garrison University, Lahore, Pakistan 

University of Engineering & technology Lahore, KSK Campus, Pakistan 

University of Education, Lahore, Vehari Campus, Pakistan 

7National Centre for Physics, Islamabad, Pakistan 


https://doi.org/10.32350/sir.73.07 


Received: November 11, 2022, Revised: February 16, 2023, Accepted: March 21, 2023, 
Published: August 28, 2023 


Hussain S, Akbar U, Ahmad M, et al. Synthesis, properties, and applications of 
carbon nanotubes: An overview. Sci Inq Rev. 2023;7(3):95-124. 
https://doi.org/10.32350/sir.73.07 


© The Authors 


This article is open access and is distributed under the terms of 
By Creative Commons Attribution 4.0 International License 


Author(s) declared no conflict of interest 


A publication of 
The School of Science 
University of Management and Technology, Lahore, Pakistan 


Synthesis, Properties, and Applications of Carbon Nanotubes: An 
Overview 


Shabbir Hussain!“, Uzma Akbar’, Muhammad Ahmad?, Muhammad Ibrar*, 
Zulfiqar Ali°, Muhammad Wagqas*, Sheikh Asrar Ahmad*°, Syed Mustansar 
Abbas’, Habib Ullah‘, and Sehrish Anwar* 


‘Institute of Chemistry, Khwaja Fareed University of Engineering and 
Information Technology, Rahim Yar Khan, Pakistan 

"Department of Chemistry, Minhaj University, Lahore, Pakistan 

3Department of Chemistry, Division of Science and Technology, University of 
Education, Lahore, Pakistan 

‘Department of Chemistry, Lahore Garrison University, Pakistan 

‘Department of Basic Sciences & Humanities, University of Engineering & 
Technology, Lahore, Pakistan 

°Department of Chemistry, University of Education, Lahore, Vehari Campus, 
Pakistan 

Nanoscience and Technology Department, National Centre for Physics, 
Islamabad, Pakistan 


ABSTRACT 


The current study attempts to review the literature concerning the 
synthesis, properties, and application of carbon nanotubes (CNTs). The 
methods used to produce carbon nanotubes include laser ablation, electric 
arc discharge, chemical vapor deposition, plasma-enhanced chemical 
vapor deposition, pulsed laser deposition, use of low-frequency ultrasound 
waves, heating a bulk polymer, and bulk sputtering. CNTs have excellent 
mechanical and thermal properties that strongly depend upon their 
structure. Functionalized magnetic CNTs are involved in magnetic force 
microscopy used in biomedicine. The liquid and plastic limit of kaolinite 
can be increased by adding CNTs to it. In the medical field, CNTs have 
numerous applications including gene delivery to cells, cancer therapy, 
drug delivery, and tissue regeneration. Their antioxidant nature also 
enables them to be used in cosmetic products and in the field of 
dermatology. They are also used to purify the environment, water, and in 
modern food-packaging technology. The sensors containing CNTs 
composite pellets are sensitive to gases, such as NH3, CO2, and CO H20. 
CNTs are used to construct gas containers for hydrogen storage. They are 
also considered ideal for structural applications and their properties can be 
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improved by making their composites with metals. Such metals may be 
introduced into the core of CNTs by different methods including solid- 
state reaction, arc-discharge method, and electrochemical techniques. The 
value of absorbed hydrogen gas in CNTs varies between 0.4 and 67 mass 
%. Recent advances encourage more research on CNTs to increase their 
clinical applications in the future. 


Keywords: carbon nanotubes (CNTs), metals, properties, synthesis 
1. INTRODUCTION 


Carbon plays an important role in almost all fields of science and 
technology, especially as a source of energy [l1]. Carbon nanotubes 
(CNTs) possess outstanding electrical and mechanical properties, which 
are used as an important component of flexible batteries [2]. They are 
hollow cylindrical structures rolled into a cylinder akin to a honeycomb 
lattice. The internal diameter of CNTs lies in the range of nanoscale or is 
lum in length. The history of carbon nanotubes (CNTs) started in 19" 
century, certainly from World War 2 when carbon fiber was first 
synthesized by Thomas A. Edison for electric light bulbs as a filament [3]. 
CNTs fall into two categories, namely single-walled carbon nanotubes 
(SWCNTs) and multi walled carbon nanotubes (MWCNTs). The 
toughness and higher thermal conductance are owed to the presence of sp” 
hybridization in CNTs. SWCNTs are insoluble in both organic and 
inorganic solvents, although they may be soluble after complexation and 
polymerization. When SWCNTs are added to aniline (organic solvent) 
then SWCNT-aniline charge transfer complex is formed that makes 
SWCNTs soluble in organic solvents and preserves their pristine nature, as 
proved by various analysis tools [4]. 


CNTs comprise graphite sheets rolled into a cylindrical pattern and 
have various applications in the fields of medicine and material sciences 
[5]. They have widespread applications in every aspect of life, such as film 
and coatings and also in the fields of biotechnology, electronics, 
environment, energy storage, and many others due to their unique higher 
thermal conductivities, toughness, and electrical conductivity [6]. 


Keeping in view the great importance of nanotechnology in various 
fields [7-10], the current study reviews the literature concerning the 
synthesis, properties, and uses of carbon nanotubes (CNTs). 
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2. SYNTHESIS OF CARBON NANOTUBES (CNTS) 

The methods used to produce CNTs include laser vaporization, electric 
arc discharge, chemical vapor deposition [11], decomposition of SiC, 
dipping graphite in cold water, torsion of graphene layers, as well as, 
mechano-thermal, pyrolysis, solar energy, liquid phase, electrolysis, and 
heat treatment of polymers [12]. 


Historically, electric arc discharge was the first method used for the 
formation of CNTs. Arc discharge method technically resembled the laser 
evaporation process. Comparison between these two methods (laser 
evaporation process and electric arc discharge) that shows that there is a 
difference in the purity and quality of the obtained products. However, the 
most favorable and developed techniques used for the synthesis of CNTs 
and related to materials on a large scale, which include different types of 
chemical vapor deposition (CVD) and arc discharge [13]. Consequentially, 
CNTs were synthesized using various techniques, such as laser ablation or 
arc discharge. However, in present days, CVD (<80°C) has replaced these 
techniques (laser ablation and arc discharge) because this method 
accurately controls the nanotube diameter, length, density, purity, 
orientation, and alignment [14]. 


2.1. Electrical Arc Discharge Method 


The arc discharge method involves the vaporization of carbon by 
applying electric field at a high-temperature gradient. This process may be 
improved in the presence of numerous metal catalysts (for instance, iron, 
cobalt, nickel, yttrium, boron, and gadolinium) under the reduced pressure 
of the inert gas [15] It also creates plasma in a glass chamber due to the 
transfer of energy from the arc to the graphite anode doped with a catalyst. 
The arc discharge setup involves two graphite rods, which act as anode 
and cathode and have a diameter of 20 mm and 7 mm, respectively. An arc 
is produced when 100-200 ampere current is provided between the 
electrodes [16]. Generally, multi-walled carbon nanotubes (MWCNTs) are 
produced when no catalysts are applied. However, single walled carbon 
nanotubes (SWCNTs) are synthesized in the presence of a transitional 
metal catalyst. The high-temperature arc discharge method (above 
1700°C) results in the formation of a mixture product and requires the 
separation of CNTs from both coal and the remaining metals (catalytic) 
[17]. The catalyst composed of nano-sized metal particles (such as Ni, Co, 
or Fe), which cause the disintegration of the precursor molecules of 
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gaseous hydrocarbon into carbon [18]. Needle-like CNTs with 
approximately 1 mm of length and 4-30 nm diameter can be produced on 
carbon cathode by using the direct arc discharge evaporation of carbon. 
Lijima et al. used the pressurized chamber filled with a gaseous mixture of 
40 torr argon and 10 torr methane [19]. The high yield production of 
CNTs depends on the concentration and nature of catalysts, gases that are 
composed of plasma, the pressure of inert gas, the arc current intensity, 
and the distance between the electrodes [17]. Figure 1 displays a 15 
kW xenon short-arc lamp used in the IMAX projection system. 


Figure 1. A 15 kW xenon Short-Arc Lamp 


Source. https://en.wikipedia.org/wiki/Arc_lamp 
2.2. Laser Ablation Technique 


CNTs can be produced using the laser ablation method (Figure 2). 
Normally, a laser is directed on the targeted carbon, which vaporizes a 
small quantity of material inside an oven warmed up to 1200°C. The 
smooth beam of laser on the target is ensured by a computer control 
system [20]. The plasma produced by this method is usually swept by 
nitrogen or argon from a high-temperature gradient and deposited onto the 
surface of the substrate, which is cooled by the external cooling system. 
Pulsed or continuous laser vaporizes 1.2% of Co/Ni with 98% composite 
of graphite target under 500 torr of the inert atmosphere of helium at 
1200°C in quartz furnace. A plume forms vapors at this high temperature 
which rapidly expands and cools. A large cluster is formed by the quick 
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compression of small carbon molecules or atoms to cool down the vapors 
[21]. The growth of nanotubes is stopped when the carbon layer cannot 
absorb particles anymore because the surface is already occupied and has 
no space for the coming particles. MWCNTs are formed by using a pure 
graphite target. The SWCNTs yield depends strongly on the metal catalyst 
type and it can be increased by increasing the temperature of the furnace, 
among other factors [12]. Laser ablation is the best method to grow 
SWCNTs with a high purity and a high quality [17]. 


Laser source 
CO», Ar-ion, 
Nd:Yag (pulsed), 


HeNe, Nd:YVO;, Sensors 
Solid state, (pyrometer, 
Diode camera, 
Raman etc.) 


(Objective) 


J 


Substrate 


(back-side laser) 


Figure 2. Laser Ablation Setup for CNT Synthesis [22] 
2.3. Chemical Vapor Deposition (CVD) Method 


In electronics and CMOS industry, chemical vapors deposition (CVD) 
is the most suitable technique for the synthesis of CNTs as compared to 
laser ablation and arc discharge methods. This technique can synthesize 
CNTs at a lower temperature [23]. CVD yields very pure and high- 
performance solid materials. This technique can be used to synthesize the 
materials in which vapor phase chemical components react at specific 
surfaces to form solid films [24]. 
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CVD was first reported in 1996 for the production of nanotubes. In 
CVD, vacuum deposition occurs for the production of CNTs. In this 
process, the substrates are exposed to a volatile precursor which reacts or 
decompose on the surface of the substrates. This method is useful for the 
production of nanotubes in a high quantity and also to control the growth 
direction onto the substrate [25]. The synthesis of CNTs through CVD 
involves the sputtering of a transition metal catalyst onto the substrate. It 
also uses thermal annealing or chemical etching to induce nucleation in 
the catalyst particle. Clusters are formed on the substrate as a result of 
thermal annealing and nanotubes are grown [26]. 


In the reaction chamber, there is a gas phase in which the carbon 
source is placed. Then, by using an energy source such as plasma, carbon 
molecule or heating coil is converted to an atomic level. Carbon 
monoxide, methane, or acetylene can be used as a carbon source [27]. 
Generally, CVD utilizes a temperature range of 650—-900°C for the 
synthesis of nanotubes, with the consequent yield of about 30% [26]. 
MWCNTs are produced by thermal CVD methods, using acetylene or 
ethylene gas as the feedstock of carbon in the presence of Ni, Co, or Fe 
nanoparticles as a catalyst. At these temperatures, carbon atoms are 
dissolved into metal nanoparticles which finally become saturated and the 
precipitation of carbon forms CNTs. The size of the metal particles 
(catalysts) can help to determine the diameter of CNTs. When other 
elements such as Mn, Cr, and Cu are used, CNTs are formed only in 
minute amounts [28]. 


The percentage yield of CNTs varies depending upon their synthetic 
methodology. Table 1 gives a brief comparison between CVD, laser 
ablation, and arc discharge methods. 


Table 1. Comparison between Different Synthesis Methods 


Arc discharge method Laser ablation method Chemical vapor deposition 
> With the help of the > The yield of the > For the synthesis of 
arc discharge tube laser ablation carbon nanotubes, CVD 
method, we can method is about is known as the most 
synthesize up to 30% 710% of SWNTs widely used method 
by weight of both with a measurable [30]. 
MWCNTs and diameter, whichis > The removal of catalyst 
SWCNTs [29]. controlled by the by acid treatment may 
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Arc discharge method Laser ablation method Chemical vapor deposition 
> The synthesis of temperature of cause destruction of the 
carbon nanotubes reaction. original structure of 
requires a > Incomparison with carbon nanotubes. 

temperature above chemical vapor > The yield may be 

1700 °C. deposition, arc increased by the use of 
> The synthesis is discharge method catalysts like MgO or 

performed in an arc is much expensive AÑO; [31]. 

discharge tube. This [29]. 


method causes less 
structural distortions 


[14]. 


2.4. Plasma-enhanced Chemical Vapor Deposition (PECVD) 


PECVD is a versatile technique that is used to obtain vertically dense 
aligned CNTs at lower temperature, as compared to the one’s used in 
CVD. It is a latest procedure for the selective positioning and vertical 
alignment of CNTs. In PECVD, the activation energy for the CVD 
reaction is provided not only by elevating temperature gradient but also by 
the energetic plasma formation in an electric field [32]. 


2.5. Pulsed Laser Deposition Method 


The CNTs can be synthesized by the pulsed laser deposition method. 
This is a thin film deposition method in which the targeted material is 
evaporated by a pulse of laser beam and a film is deposited on the 
substrate surface. The furnace contains the targeted substance placed at the 
bottom in which the substrate is placed at the top side of the furnace. The 
laser beam normally used is Nd: YAG. The laser beam strikes the targeted 
atom and vaporize them. The vaporized atoms are called a plume. The 
plume is moved towards the substrate, deposited and grows as carbon 
nanotubes upon the substrate surface [33]. The rate of deposition and the 
laser beam can be collaborated to control certain parameters of the CNTs. 
These CNTs are deposited in the substrate surface that can be separated 
from the carbon nanoparticles, amorphous carbon, and other unwanted 
impurities. These prepared CNTs can be purified by several methods that 
includes gas-phase purification and liquid phase purification, which are 
used most commonly for the purification procedure [34]. 
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2.5.1. Gas-phase Purification. In the gas-phase purification method, 
CNTs are purified by applying high-temperature oxidation, which is 
continued by the repeated extractions with nitric acid and HCl. In this 
process, the synthesized CNTs have a high purity and stability, including 
little impurities or residual catalyst [34]. 


2.5.2. Liquid phase Purification. In liquid phase purification, the 
CNTs are purified in several steps, which mainly include: 


i. The preliminary filtration, which is done to remove the bulk residual 
graphite particles. 


ii. Then the CNTs are dissolved in both conc., acid, and organic solvents 
to remove the unreacted catalyst and fullerenes. 


iii. Next is the centrifugal separation of CNTs from the solution. 
iv. Then microfiltration is done. 


v. The final step of purification is chromatography, which is used to 
separate the MWCNTs and SWCNTs [34]. 


2.6. Modification of CNTs by Low-Frequency Ultrasound Waves 


Multi-walled carbon nanotubes (CNTs) can be modified by treating 
with 20 kHz ultrasound in combination with dilute HNO3 and dilute 
H2SO0; for 30 min at 12 W cm” sonication. These specified conditions 
prevent the aggregation of the nanotubes and allow an efficient dispersion 
in ethanol or in chitosan [35]. 


2.7. Synthesis of CNTs through Bulk Polymer 


The CNTs can be synthesized by heating a bulk polymer at about 400 
°C in the air. This can be obtained by polyesterfication reaction of 
ethylene glycol and citric acid. The formation of CNTs is confirmed by 
the number of spectroscopic analyses. The lengths of CNTs is generally 
less than 1 um, whereas their diameter ranges from 5-20 nm [36]. 


2.8. Synthesis of CNTs through Bulk Sputtering 


One of the most promising methods among all is bulk sputtering 
method, which may be employed for the mass production of CNTs, 
however, this process is overpriced[37]. During the process of sputtering, 
the targeted graphene is usually placed as the targeted material in the 
chamber. The reaction chamber is filled with an inert argon gas, which is 
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commonly used when an inert atmosphere is required or for the production 
of titanium and other reactive material. The electric field is applied 
between the two electrodes, which accelerate the ionization process of 
argon gas to produce the charged particles that strike the surface of the 
targeted graphene material and eject the carbon atoms from the surface of 
graphene [38]. The CNTs are grown on the surface of the substrate. Thus, 
the CNTs prepared by this method are impure and can be purified by 
several other methods [37]. 


3. PROPERTIES 


Carbon nanotubes (CNTs) have unique properties due to which they 
are found in numerous applications, which have a great impact on human 
health. However, some chemicals have also evolved during CNTs 
production and handling, which are usually inhaled during the inhalation 
process [39]. The tensile strength of CNTs is hundred times greater than 
their steel and thermal/electrical conductivities, which are comparable to 
copper [25]. CNTs show significant properties including physicochemical 
properties, metallic or semi-metallic characteristics, great mechanical and 
electrical characteristics, good thermal conductivity, high electrical 
conductivity, and a large surface area [40]. They are considered ideal for 
the production of next-generation composite materials due to their specific 
properties [41]. The properties of CNTs can be modified by doping with 
heteroatoms such as nitrogen and boron to tune their physicochemical 
characteristics for specific applications [42]. 


3.1. Electronic Properties 


CNTs possess hollow cylinders with approximately 20 carbon atoms 
around the circumference of the cylinders and microns in length. They can 
act as semiconductors or conductors depending on their structure and 
possess characteristic electronic properties. CNTs have unusual electronic 
properties due to their high electrical conductivities as compared to 
copper. The symmetry of the planner system can be broken down by 
rolling action, which causes a specific direction with respect to the 
hexagonal lattice and axial directions [43]. Different theoretical and 
experimental works suggest that single-walled CNTs are considered as 
quantum wires, which have one dimension [44]. The graphene sheet is the 
main site for single-walled CNTs having a main electronic structure and is 
considered to possess metal properties with bands of conduction [45]. 
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Moreover, carbon nanotubes are also considered as an example of ideal 
quantized one-dimensional conductors [46]. 


3.2. Mechanical Properties of CNTs 


CNTs are one-dimensional materials that display the mechanical 
properties, such as the tensile strength and young’s modulus. However, it 
is difficult to produce the pure form of CNTs and to control and manage 
their properties. Due to these properties, CNTs are considered good for 
polymer enforcement. The extensive mechanical characteristics of CNTs 
are due to the sp? hybridization of carbon-carbon double bond. The 
densities of carbon nanoparticles are very low (1.3 gem™) as compared to 
stainless steel. Young’s modulus of CNTs is higher than carbon fibers 
with larger values than that of 1TPa, which is about five times larger than 
steel [43]. CNTs have excellent mechanical properties that strongly 
depend upon their structure. There are three kinds of nanotubes (SWNTs, 
MWNTs, and catalytic MWNTs) studied under laboratory actions that 
always consist of structural defects [47]. The calculation for the stiffness 
constant of SWNTs can be done by adopting the elastic modulus of 
graphite to understand the important properties like mechanical properties 
of graphite of a single crystal [48]. 


3.3. Thermal Properties of CNTs 


One of the most important properties of CNTs is thermal property, 
which is directly related to their small sizes and unique structures. 
Thermal and specific heat conductivity of bulk MWNTs can help to 
measure the thermal properties of carbon nanoparticles [49]. The thermal 
characteristics of CNTs are directly related to each other’s as the size of 
devices (electronic or mechanical) are reduced to the micro and nanometer 
level. With the help of molecular modeling predictions, the studies of 
nano-electronic devices became very easy, however, it is quite difficult to 
measure the accurate thermal properties [50]. The excellent thermal 
conductivity of 3000 W/m K at room temperature of SWNTs makes them 
important for thermal applications [51]. The high thermal conductivity of 
nanotubes has a number of applications in thermal management such as 
heat sinking of silicon processors, or increasing the thermal conductivity 
of plastics in areas of housing for electric appliances [52]. Carbon 
nanotubes show greater variations in thermal conductivities and the 
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conductivity of nanotubes is greater along the axis due to greater 
anisotropic property [53]. 


3.4. Magnetic Characteristics and Heat Generation 


CNTs are widely accepted due to their fundamental properties. The 
basic field interest of CNTs is the magnetic functionalization, which has a 
wide applications as nano-scaled magnetic objects [54]. The 
semiconducting and metallic properties of CNTs can be predicted by both 
first principal calculation and tight binding that depends upon their size 
and helicity. Few unusual characteristics were demonstrated by many 
recent experiments [55]. For a magnetic field, high diamagnetic sensitivity 
is found in both perpendicular and parallel axis of the tube [56]. 
Specifically, the magnetic functionalization gives tremendous potential to 
CNTs that can provide a practical approach to give a stable coating for 
protection from degradation and oxidation. The functionalized magnetic 
CNTs are involved in magnetic force microscopy, which is used in 
biomedicine in the form of magnetic nano-vectors or spintronics [57]. 


3.5. Carbon Nanotubes (CNTs) Composites 


Due to the magnificent physical properties of CNTs, they form a 
metal matrix with the metals to take the advantage of their high tensile 
strength and electric conductivity. They are also considered ideal for 
structural applications and their properties can be improved by making 
their composites with metals. The following techniques are used for the 
fabrication of CNTs polymer nanocomposite material. 


3.5.1. Solvent Casting. It is an important procedure for the preparation 
of CNTs concentrated polymer nanocomposite materials. These solvents 
have a significant influence on the properties of nanocomposites. The in 
solvent casting of the nanotubes is facilitated by dispersion and it involves 
the preparation of a suspension of CNTs in the solution of a desirable 
polymer by energetic agitation. The process involves the dissolution of a 
polymer in an organic solvent and then the addition of particles (generally 
salts) with specific dimensions to the solution. It is followed by the 
shaping of consequent mixture, its final geometry a membrane that can be 
produced by casting it onto a glass plate or a scaffold is formed by using a 
three-dimensional mold. After evaporation of solvent evaporates, a 
composite material is produced, which contains the particles together with 
the polymer [58]. 
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3.5.2. Melt Mixing Method. CNTs can also be synthesized by the 
melt mixing method, which is used for the thermoplastic polymers that 
gets soften on heating. This process involves an elevated temperature 
gradient to decrease the viscosity of the substrate and high shear forces to 
the nanotubes bundle. Templates of numerous shapes can then be obtained 
by various fabrication techniques, including injection molding, 
compression molding or extrusion formed composites of commercial 
polymers, such as acrylonitrile butadiene styrene (ABS) with MWCNT, 
polypropylene, and high impact polystyrene [59]. 


For the achievement of multifunction and high-performance stability, 
CNTs are considered to be an ideal filler for polymer type matrixes, due to 
their sizes in nanometers, high and well-known aspect ratio, and more 
specifically due to their amazing strength and high thermal and electrical 
conductivities [60]. Solution casting and melt blending are the most 
common methods that are used to produce polymer composites; inorganic 
fillers may also be used in making carbon nanotubes-polymer composites 
[61]. In recent times, fabrication of CNTs with the polymeric substances 
gives the characteristics synergistic effects to the host and boosts a vast 
interest in the development of multifunctional CNTs along with alternative 
materials [62]. 


4. APPLICATIONS OF CARBON NANOTUBES (CNTS) 


Nanoparticles are found in a broad range of biological [63—66] and 
non-biological applications [67-70]. CNTs functional materials 
(CNTFMs) had shown a potential impact on the fields of science, 
engineering, and technology with the transformation of nanoscience 
towards a practical applications [71]. 


4.1. Carbon Nanotubes as Stabilizers 


The liquid and plastic limit of kaolinite can be increased by the 
significant addition of CNTs to it. The resultant mixture would have lower 
soil strength, higher compressibility, and reduced hydraulic conductivity 
[72]. The presence of CNTs in construction materials increase the 
mechanical strength of a material. It has the capacity to increase the 
flexibility of a material [73]. It also decrease the breakability of a material. 
Soil cement itself is not a flexible material and it can be easily broken by 
applying stress. Consequently, materials which are made by only soil- 
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cement can be easily breakable by applying any kind of stress. CNTs have 
significantly solved such type of problems [74]. 


4.2. CNTs Applications in the Medical Field 


In the medical field, CNTs have numerous applications including gene 
delivery to cells, cancer therapy, drug delivery, and tissue regeneration. 
MWNTs and SWNTs have a greater potential to improve the traditional 
drug delivery to the cells [75]. CNTs as compared to other nanocarriers 
that can be easily modified for conjugation of bioactive compounds and 
ligands for the targeting process. Additionally, they also found diagnostic 
applications to drug delivery [76]. The use of CNTs in the nerve tissue are 
explained in several reports [77]. 


CNTs applications are also found in drug delivery, genes, cells, and 
cytokines. Moreover, they can also develop tissue induction and cell 
activation properties along with futuristic biomaterials characteristics, 
which are beneficial for the human body. For the self-repaired potential of 
the human body, it is very important to use CNTs flexibly to adapt to their 
environment [78]. Their solubility is applicable in biocompatibility, 
secretion, blood transportation, and gastrointestinal absorption. The CNTs 
composites are also involved in therapeutic drug delivery systems [79]. 
Significantly, CNT dispersions should have a uniform and stable 
distribution to a sufficient degree for an accurate concentration. The 
solubility of pristine CNTs in the aqueous solution is an important factor 
for their use as a practical drug carrier, which can be owed to the 
hydrophobic property of the graphene sidewalls and also m-n strong 
interaction between CNTs. There had been a prime focus on nano-carbon 
biomaterial research and its development over the last 15 years. However, 
it is strongly believed that this focus would lead to paradigm shifts and 
lead to major advancements in global medicine by following the road 
towards clinical applications [80]. 


Figures 3 and 4 display important applications of CNTs in the field of 
medical. 
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Figure 4. Application of CNT in Biomedical Field [82] 
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4.3. Application of Carbon Nanotubes (CNTs) as a carrier of the 
Anticancer Drug 


Chemotherapy is usually coupled with other treatment techniques to 
reduce the size of the tumor in traditional methods, however, this 
technique also causes toxicity to the other living cells, which has adverse 
side effects. To reduce the side effects, a new method has been developed 
in which anticancer drug is delivered to tumors by CNTs [83]. To deliver 
the drug to the target cell, CNTs are functionalized with the chemical 
receptor (protein or nucleic acid), and then the anticancer drug is fixed in 
functionalized CNTs with the open ends; this drug along with CNTs is 
introduced into the animal body either by oral or injection. Finally, drug 
carrier CNTs pass across the nuclear membrane and cytoplasmic 
membrane without generating toxic effects and reach the targeted cell. 
There are two targeted ways through which a drug is delivered to a cell; 
one is non-internalization, while the second one is the internalization 
method. Comparatively, the second method is considered better than the 
first method. In the internalization method, both drug and CNTs are 
entered inside the cell from the intracellular environment that assists to 
deteriorate the drug carrier conjugate for releasing drugs into the cell. 
While in the method of non-internalization, the extracellular status helps 
to deteriorate the drug carrier to conjugate, then the drug itself penetrates 
across the lipid membrane to enter into the cell; here the drug is degraded 
before reaching the targeted cell. There are two possible mechanisms of 
internalization of CNTs, which include insertion and diffusion and 
endocytosis. The whole procedure is used for cancer therapy [84]. 


4.4. Carbon Nanotubes as Preservatives 


Many investigations have focused on the development of novel 
antioxidants due to their importance in food and pharmaceutical industry 
[85]. CNTs and nano-horns are antioxidants in nature. They can be used in 
drugs to reduce the process of oxidation in the human body. Their 
antioxidant nature also enables them to be utilized in cosmetic products 
and in the field of dermatology along with zinc oxide sunscreen, which is 
used to avoid oxidation of vital skin constituents [86]. For convenient and 
processed food, aromatic organic acids are used to inhibit a vast range of 
fungi, molds, bacteria, and yeasts [87]. Despite this, extensive 
consumption of these types of preservatives in food may be harmful for 
human health and can cause several allergic diseases like dermatitis, hives, 
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convulsion, and others in humans [88]. Currently, carbon nanomaterials 
have attracted significant attention worldwide in modern food-packaging 
technology due to their easy functionalization, high surface area-to- 
volume ratio, excellent antimicrobial activity, superb physico-mechanical 
and water resistance properties, strong adsorptive ability, high thermal, 
and electrical conductivity,, CNTs-based nano-sensors offer a better 
assessment of the freshness, security, safety, quality, and packaging of 
food products [89]. 


4.5. Application of Carbon Nanotubes (CNTs) as a Sensor 


Sensors are mostly imported as detecting devices, which are used in 
different fields. There are different types of sensors such as biosensors and 
molecular sensors. The efficiency of sensors can be enhanced by attaching 
CNTs to them (Figure 5) [90]. CNTs have been employed in the sensing 
and detection of liver and pancreas cancer [91]. The usages of CNTs for 
the detection and sensing of NO2 gas are well recognized [92]. CNTs 
incorporated sensors can be used to bring revolutionary changes in 
different sectors, especially in the biomedical industry [83]. The sensors 
containing CNTs composite pellets may be constructed by attaching 
different chemical groups onto the end; they are also sensitive to gases 
such as NH3, CO2, and CO H2O [93]. Gas sensor are used in medical, 
industrial, and in commercial areas for the reduction of greenhouse gases 
and to monitor the environment of the combustion engine [94]. Ammonia, 
carbon dioxide, and oxygen gas can be detected using their conductivity 
and permittivity properties by multi-walled nanotubes [95]. 


The CNTs have unique physical, adsorption, and electrochemical 
properties. The strong adsorption capacity of the CNTs and their good 
sensitivity towards atoms and molecules adsorbed on their surface enable 
them to design sensors-based CNTs. Many gas sensors (detectors) based 
on the CNT have been reported in previous literature [96] that include: 


i. Shift gas sensors 

ii. Resonance frequency 
ili. Capacitance gas sensors 
iv. Ionization gas sensors 


v. Adsorption gas sensors 
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Their main operating principle involves adsorption during which an 
adsorbed gas molecule transfers an electron to or takes it from a nanotube. 
This changes the electrical properties of the CNTs, which can be detected 
and measured. There are numerous gas sensors based on pure SWCNTs 
and MWCNTs modified by metals, metal oxides, polymers or various 
functional groups [97]. 


PRESSURE 


Figure 5. Application of CNT as a sensor [98] 


4.6. Application of Carbon Nanotubes (CNTs) in Water Filtration and 
Environment 


Carbon nanotubes (CNTs) find numerous applications in the purification 
of water [99] and the remediation of pollutants [100]. They can reduce the 
distillation cost by up to 75% and their membranes help in filtration. 
CNTs have a very small particle size, which enable smaller particles (like 
water molecules) to pass through them, while larger particles are blocked 
(they cannot pass through carbon nanotubes), such as chloride salts [101]. 


4.7. Carbon Nanotubes (CNTs) in Electronic Devices 


In electronic devices, CNTs can be used as field emission sources, 
which is done by applying the potential between CNTs surface and anode. 
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Due to the curvature present in the form of pentagons in CNTs , electrons 
are discharged from their tips [102]. The use of CNTs as electron emitters 
is associated with several advantages, which include stable field emission, 
high current densities over prolonged periods, low emission threshold 
potential, and long lifetime of the components for the construction of field 
emission devices [103]. 


CNTs can be used in supercapacitors because they have a large 
surface area [104]. Therefore, by using carbon nanotubes (CNTs), the use 
of platinum can be reduced, which creates some problems specially in 
fuel cells [105]. 


4.8. Gas and Hydrogen Storage 


Carbon nanotubes (CNTs) are used as a metal container because of 
their hollow cylindrical nature. Metals may be introduced into CNTs core 
by different methods, including solid-state reaction, arc-discharge method, 
and electrochemical techniques [106]. The environmental problems are 
rising day by day due to the overuse of fossil fuels. There are various 
primary energy sources, such as, wind, solar energy, thermonuclear 
energy, and geothermal. The best choice among all is hydrogen because 
of its adverse advantages. It can be easily produced, it has a high 
utilization efficiency and it transforms without any additional problems at 
the consumer end [107]. CNTs are used to construct gas containers, which 
are used to store hydrogen gas. Hydrogen gas is absorbed in single-walled 
nanotubes. The value of absorbed hydrogen gas in CNTs varies between 
0.4 and 67 mass % [82]. 


5. CONCLUSION 


Carbon nanotubes can be produced by laser ablation, electric arc 
discharge, chemical vapor deposition (CVD), plasma-enhanced chemical 
vapor deposition (PECVD), pulsed laser deposition, use of low-frequency 
ultrasound waves, by heating a bulk polymer and bulk sputtering etc. 
Chemical vapor deposition (CVD) and arc discharge are the most favorite 
methods for CNTs production on large scale; however, the latter is 
comparatively expensive. The CVD yields very pure and high- 
performance solid materials and has the ability to control accurately the 
nanotube diameter, length, density, purity, orientation, and alignment. 
Plasma-enhanced chemical vapor deposition (PECVD) is a latest 
technique for the selective positioning and vertical alignment of CNTs. 
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The laser ablation is the best suitable method to grow single-wall 
nanotubes with high purity and high quality. The CNTs prepared by 
pulsed laser deposition require either gas- or liquid-phase purification. 
CNTs show characteristic metallic/semi-metallic, mechanical, electrical, 
thermal and electrical properties which are improved by making their 
composites with metals. They can improve the mechanical strength and 
flexibility of a construction material. CNTs find applications in diagnosis, 
drug delivery, gene delivery to cells, cancer therapy, drug delivery, tissue 
regeneration, biomedicine, cosmetic products, dermatology, 
environmental science, water purification, modern food-packaging 
technology, electronic devices, gas sensors and in construction of gas 
containers for hydrogen storage. 
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